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ABSTRACT 

We study the physical and photometric properties of galaxies ed z — A in cosmological hydrodynamic 
simulations of a ACDM universe. We focus on galaxies satisfying the "B-dropout" criteria of the 
Great Observatories Origins Survey (GOODS). Our goals are: (1) to study the nature of high-redshift 
galaxies; (2) to test the simulations against published measurements of high-redshift galaxies; (3) to 
find relations between photometric measurements by HST/ ACS (0.4 - 1 fim) and Spitzer/IRAC (3.6 
- 8 /im) and the intrinsic physical properties of GOODS "B-dropouts" such as stellar mass, stellar 
age, dust reddening, and star-formation rate; and (4) to assess how representative the GOODS survey 
is at this epoch. Our simulations predict that high-redshift galaxies show strong correlations in star 
formation rate versus stellar mass, and weaker correlations versus environment and age, such that 
GOODS galaxies are predicted to be the most massive, most rapidly star-forming galaxies at that 
epoch, living preferentially in dense regions. The simulated rest-frame UV luminosity function (LF) 
and integrated luminosity density are in broad agreement with observations at z ~ 4. The predicted 
faint end slope is intrinsically steep, but becomes shallower and is in reasonable agreement with data 
once GOODS selection criteria are imposed. The predicted rest-frame optical (observed 3.6^) LF 
is similar to the rest-frame UV function, shifted roughly one magnitude (AB) brighter. We predict 
that GOODS detects less than 50% of the total stellar mass density formed in galaxies more massive 
than log(M*/M0) > 8.7 by z = 4, mainly because of brightness limits in the HST/ ACS bands. 
Most of these results are somewhat sensitive to the prescription used to model the effects of dust 
extinction. We develop a physically-motivated model that is based on taking the simulation-predicted 
metallicities and using the dust-metallicity relation calibrated locally from SDSS. This model generally 
produces results in good agreement with observations, although it produces a modest excess of bright, 
rapidly star-forming galaxies. The slope of the predicted stellar mass-metallicity relation is in excellent 
agreement with low-redshift measurements of the stellar mass-gas-phase metallicity relation such that 
over two decades in stellar mass, galaxies are less enriched than low-redshift galaxies of similar stellar 
mass by roughly 0.6 dex. The most rapidly star forming galaxies in our simulations have rates 
exceeding 1000 Mq yT~^, similar to observed sub-mm galaxies. These galaxies are not starbursts, 
however, as their star formation rates show at most a mild excess (~ 2 — 3x) over the star formation 
rates that would be expected for their stellar mass. It is possible that the observable counterparts to 
these bright galaxies do not follow our dust prescription and are instead heavily extinguished. The 
overall distribution of dust reddening and mean stellar age may be constrained from color-color plots 
although the specific value for each galaxy cannot. 

Subject headings: cosmology: theory — galaxies: evolution — galaxies: formation — galaxies: high- 
redshift — galaxies: photometry — galaxies: stellar content 



1. INTRODUCTION 

According to the currently-favored hierarchical model 
of structure formation, the first generations of stars be- 
gan forming in low-mass galaxies at very high redshift. 
These galaxies then grew by accreting gas from their sur- 
rounding intergalactic medium (IGM), and merged into 
larger galaxies that we observe today. The details of 
how this process operates are still not well understood. 
Numerical simulations that include a range of physical 
processes believed to govern galaxy formation allow the 
hierarchical scenario to be tested through detailed com- 
parisons with observed galaxies. This yields insights into 
the nature of galaxy formati on and assists wit h interpret- 
ing observations of galaxies ijWeinberg et alJ l20G2). 

^ University of Arizona, Department of Astronomy, Tuscon, AZ 
85721 

^ Harvard-Smithsonian Center for Astrophysics, 60 Garden 
Street, Cambridge, MA 02138 



The current generation of deep surveys from the opti- 
cal to the radio are giving us an unprecedented view of 
the high-redshift {z > 3) universe, and as a result are 
providing new and stringent tests on models of galaxy 
formation. Recent years have also seen a rapid growth in 
the power and sophistication of numerical simulations of 
galaxy formation, so that detailed comparisons are ever 
more meaningful and illustrative. For these reasons, the 
time is ripe for a detailed investigation into the physical 
nature and observable properties of high-redshift galaxies 
in cosmological simulations of galaxy formation. 

The most effective method for finding high-redshift ob- 
jects for o bservational study is the Ly man break selection 
technique ijSteidel fc Hamilto n 1993; Steidel et al. 1996), 
exploiting the large differential in flux across the rest- 
frame Lyman limit owing to internal absorption within 
the galaxy as well as strong absorption by intergalactic 
neutral hydrogen along the line of sight. Together with 
the fact that high-redshift galaxies are necessarily young 
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and therefore blue, this endows high-redshift galaxies 
with distinctive colors that make photometric selection 
straightforward. Spectroscopic foUowup has shown that 
this technique is highly effective at isolating Lyman break 
galaxies (LBGs) while suffering little contamination from 
lower redshifts. 

The first attempts to simu late LBGs used semi- 
analytic models (SAM s; see e.g. iKauffmann et al.lll993t 
ISomerville et aljl2001() . which employ merger trees or 
gravity-only numerical simulations to follow the growth 
of dark matter perturbations, and use analytic prescrip- 
tions to account for the remaining physical processes 
involved in forming galaxies such as merging, cooling, 
and feedback. These prescriptions give SAMs flexibil- 
ity and low computational cost, but their predictions of 
galaxy properties consequently rely on a variety of sim- 
plifying assumptions. Hydrodynamic simulations com- 
plement the SAM approach by allowing galaxies to be 
simulated with minimal assumptions regarding the dy- 
namics of gas and dark matter, but are computationally 
costly and are limited by numerical resolution and vol- 
ume effects. 

Different SAM prescriptions have yielded different 
physical models for L BGs, with some SAMs (notably 
ISomerville et aljl2001|) indicating that a significant frac- 
tion o£JLBGs_are merger-induced starbursts, while oth- 
ers (jBaugh et al.lll99^ holding that LBGs are the most 
massive galaxies to have formed at these epochs. Hy- 
drodynamic simulations, meanwhile, have consistently 
Dointed towards the latter scenario fDave et al." 1999; 
iNagamine ct al. 2004a, 20G5a; Weinberg et al. 2002). To 
date, efforts to infer the physical properties of LBGs via 
population synthesis modeling have tended t o support 
the massive galaxy hypothesis ( Shap lev et al.l2 0G5. 2001; 
iPapovich etaI]l200ll:lBarmbv et alil2004H although only 
a small subset of LBGs have had their physical properties 
studied in detail. In this paper, we elaborate on the pre- 
dictions of the massive galaxy model in order to facilitate 
further comparison with current and future observational 
work. 

While determining the basic properties of LBGs as 
a constraint on galaxy evolution models is an impor- 
tant goal unto itself, it is also a critical step towards 
constraining specific volume-averaged quantities such as 
galaxy number densities, luminosity functions (LFs), lu- 
minosity d ensity (pi/), a nd the history of cosmic star 
formation (|Madau et al.1 119961: ILillv 1 11997^ . All mod- 
els of galaxy formation predict a substantial popula- 
tion of sources at high redshift that cannot be detected 
via the Lyman break technique, although they differ 
as to t his populat i on's p roperties. For example, the 
SAM of lldzi et all (20M, Figure 2) predicts a correla- 
tion between rest-frame UV flux and stellar mass that 
is substantially we aker than t h e tren d from the hydro- 
dynamic model of IDave et al.l ljl999l Figure 3; see also 
our F igure [Hjl . Thus, the SAM used by Id zi et all 
(|200^ predicts that UV-faint and UV-bright galaxies 
are more similar in stellar mass — implying that a larger 
fraction of total formed stellar mass may be "missed" — 
th an would be pred icted by the hydrodynamic model 
of IDave et alJ ljl999(l . Testing the observable properties 
of LBGs against models will help in accurately assessing 
their contribution to global stellar mass and star forma- 
tion rate budget, as well as their relationship to various 



other galaxy populations at those epochs. 

Currently the single greatest difficulty in modeling 
LBGs accurately is the poorly-understood impact of ex- 
tinction owing to dust. Nearly all star- forming galax- 
ies at high re dshift are seen to suffer fr om some level 
of ex tinction (IjAdelberger fc Steidel 200fli: FShaplev et alJ 
12001(1 . usually expressed as a color excess E{B — V). To 
date, constraints on reddening at high redshift have been 
obtained by fit ting population synthesis models to ob - 
served LBGs (Shap lev et al.ll200ltlPapovich et alJl200H) . 
Unfortunately, this means the resulting distributions of 
E{B — V) are sensitive to selection biases caused by the 
reddening itself, making it difficult to infer how many 
galaxies may be missed by such studies. In particu- 
lar, as LBG selection criteria are biased against heav- 
ily re ddened objects with E{B - V) > 0.3 (Dad di et alJ 
|2004|) . LBG studies are likely to underpredict the me- 
dian E{B — V) suffered by the complete galaxy pop- 
ulation. In support of this, recent observations sug- 
gest that a substantial fraction of the total star forma- 
tion occurring at high redshift occurs in highly reddened 
galaxies that are too faint in rest-frame UV to be se- 
lected as LBGs but c an be selected eithe r via their rest- 
frame optic al colors l|Dadd i et a ^ |2004[) or in the sub- 
millimeter (|Chapman et al.ll2004l |2005|) . The effects of 
dust on the most massive galaxies (log(M»/M0) > 11) 
m ay not even follow simple attenuation laws su ch as that 
of lGalzetti et al.l l|2000|) l|Chapman et alJl2005l: Papovich 
et al. 2005, in preparation). Insight into the effects of 
dust on a galaxy's observed spectral energy distribu- 
tion (SED) and the relationship of extinction to its more 
fundamental properties such as stellar mass, metallic- 
ity, and star formation rate is necessary for determining 
the completeness of the current observational census of 
high-redshift galaxies. In this paper we study a range 
of reddening models, comparing to observations where 
possible, and provide an observationally-motivated red- 
dening prescription that is the most sophisticated one 
used to-date in simulations. 

In order to perform a meaningful statistical compari- 
son to understand these issues, a large sample of observed 
LBGs is required. The Great Observatories Origins Deep 
Survey (GOODS) has provide d an unprecedented view 
into the high-redshift universe ((Giavahsco et alJf2004a^ . 
GOODS probes to a depth {iab < 27.5) comparable 
to previous space-based observat ions such as the Hubble 
Deep Field ((Winiams et alJll996D while covering an area 
(320 arcmin^) comparable to shallower ground-based 
studies. Combining HST/ ACS data and Spitzer/IRAC 
data allows constraints to be placed on the star forma- 
tion rate, stellar mass, extinction, and other properties 
of galaxie s out to z ^ 4 an d beyond l|Giavalisco et alj 
l2"004k .Papovich et al.ll2004|) The field size is e ven large 
enough to study clustering l)R,oche et al.1 1200.'!) . and X- 
ray observations have provided aux iliary da ta on AGN 
activity and star formation (Lehm er et alJ |2004) . At 
present, 1115 LBGs at z ^ 4 have been isolated from 
GOODS data ( Giavali sco et al.ll2004b) . Our simulations 
have sufficient dynamic range to numerically resolve the 
entire "B-dropout" sample within a comoving cosmolog- 
ical volume that is comparable to the volume probed 
by GOODS at z ~ 4. In addition, these simulations 
have been shown to broadl y match the observed c osmic 
star formation rate density l|Springel fc Hernouistl 
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iHernauist fc Springe]|l200^ . and the rest-frame UV lu- 
minosity functions of LBGs when a moderate amount 
of dust ex tinction is ass umed ijNaeamine et al.ir2004albl 
|l)05a; Ni ght et al.ll20nl . 

In this work we present comparisons between the en- 
semble statistical properties of high-redshift galaxies in 
our simulations and in the GOODS data set; examine the 
basic physical properties of a simulated GOODS sample 
and discuss the trends that emerge; investigate how the 
photometric properties of the GOODS sample relate to 
the underlying physical properties; determine how the 
model parameters affect the observed color-magnitude 
and color-color relations; discuss the extent to which 
the observed GOODS sample will be representative of 
the complete galaxy population at z = 4; and compare 
our results with those of previous simulations and semi- 
analytic models. Our work provides a thorough investi- 
gation into using the observed photometric properties of 
high-redshift galaxies as a new frontier for testing galaxy 
formation models. 

The paper is organized as follows: In Section |21 we 
summarize the input physics in the simulation, our pro- 
cedure for identifying galaxies and computing their ob- 
served photometry, and our procedure for adding dust 
reddening to the simulated galaxies. In Section we 
discuss several tests that we used to verify that our sim- 
ulated galaxies' properties are numerically resolved. In 
Section ^ we summarize the ensemble statistical prop- 
erties of the simulated sample including incompleteness 
owing to the GOODS z ~ 4 selection function, the in- 
tegrated luminosity density, and LFs. In Section [S] we 
present a number of plots comparing physical and photo- 
metric properties and briefly introduce the most rapidly 
star forming galaxies in the simulated sample. Finally, 
in Section we present our conclusions. 

2. SIMULATIONS AND SAMPLE DEFINITION 

2.1. Simulations 

Our work is based on cosmological hydrodynamic s 
simulations done using GADGET-2 l|Springe]l l200fiD . 
GADGET-2^ models the evolution of a system of parti- 
cles under gravity using a Tree-particle-mesh solver and 
hydrodynamical forces using an entropy-conservative for- 
mulation ( Springcl & Hcrnquist 2002) of smoothed par- 
ticle hydrodynamics (SPH). Our particular simulation is 
the G6 si mulation, an extension of the G -series runs de- 
scribed in lSnringel fc HernouistI l|2008bD . It was started 
at z = 79 with 486^ dark matter and 486^ SPH par- 
ticles in a cube with periodic boundary conditions and 
comoving side length lOO/i^^ Mpc. The simulation has 
a SPH particle mass resolution of 9.7 x 10^ h''^ Mq 
and an equivalent plummer softening length of 5.33 
kpc. For comparison, we also employ the D5 simula- 
tion of Soringcl & Hcrnquist (2003b), which has higher 
resolution than G6, but in a smaller volume (see be- 
low) .Tl]B^;_wer£mn_using a concordant ACDM cosmol- 
ogy llSnergel et a,lJl2nnl: (^ini, ^a, as, h) = (0.3, 
0.7, 0.04, 0.9, 0.7); we assume this cosmology throughout. 
These simulations were run on the Athlon-MP cluster at 
the Center for Parallel Astrophysical Computing (CPAC) 
at the Harvard-Smithsonian Center for Astrophysics. 

^ Pu blically available at ^^^^^^^ 
http : //www . mpa-garching . mpg . de/gadget/| 



Galaxy star formation histories are the key output 
of GADGET-2 that we use to model the observable 
properties of high redshift galaxies. GADGET-2 forms 
stars using a subgri d multiphase model for the inter- 
stellar medium flSM: ISnringel & Hernaui st' 2003a). If a 
SPH particle's density surpasses a threshold value based 
on the local thermal Jeans mass, it acquires a two- 
phase ISM consisting of a hot ambient medium in pres- 
sure equilibrium with a cold dense medium (based on 
iMcKee fc Ostrikedll977ll . A thermal instabihty causes 
gas in the hot phase to cool and condense into the cold 
medium. The cold medium, in turn, forms stars on a 
timescale that is calibrated to match the obse rved star 
formation rate in local spirals l|Kennicutt)l 1 998|) . A frac- 
tion of these stars forms supernovae that enrich and de- 
posit thermal energy into the hot phase of the ISM. The 
thermal feedback causes cold gas to evaporate back into 
the ambient hot medium, depleting the gas reservoir from 
which stars form. In this way, GADGET-2 computes a 
subresolution, self-regulated star formation and feedback 
cycle within each star-forming SPH particle. The code 
also tracks the metal enrichment of each SPH particle 
as the simulation evolves. Star particles are spawned in 
two steps from each SPH particle by GADGET-2, so that 
each star particle has half the mass of the SPH particles. 
A star particle inherits the metallicity of its parent SPH 
particle, and does not change thereafter. 

Unfortunately, by itself this multiphase ISM model 
cannot prevent the simulation from overproducing the 
observed stellar mass density in the imi vers e, so 
some form of kinetic feedback is needed (Balog h et alJ 
2m[ IKav et all 12 002: Soring el fc HernauistI I2003at 
Hernauist fc Snringel .2003.) . For this reason, the sim- 
ulation also includes a phenomenological Monte Carlo 
prescription for generating galactic superwinds as ob- 
served in local starburst galaxies. The code selects SPH 
particles with a star-forming ISM and gives a "kick" 
to their velocity vectors in the direction of the cross 
product of the particle's velocity and acceleration, cre- 
ating broadly bipolar outflows. The model assumes 
that the mass outflow rate from star forming regions 
is twice the local star formation rate, and the veloc- 
ity of each "kick" is set to 484 km s'-^ in order to re- 
produce the observed space density of s tellar mass at 
low redshift ijSpringel fc Hernauistil2003b|) . Besides sup- 
pressing star formation, this model creates qualitatively 
realistic large-scale galactic outflows that chemically en- 
rich th e IGM as proposed bylAguirre et ajj t,2001a.bil , al- 
though lAguirre et al.l l)2005(l suggest that this superwind 
model may require some tuning. 

2.2. Galaxy Identification 

We use the Spline Kernel Interpolative DENMAX 
(SKID)'' group finder to identify gravitationally bound 
groups of star and gas particles as galaxies. SKID oper- 
ates in several steps, which may be summarized as fol- 
lows: (1) slide particles along the gradient of the initial 
baryonic density field until they are confined into sub- 
groups about the initial density peaks; (2) combine sub- 
groups into groups using a friends-of-friends algorithm; 
(3) reject particles from each group that are not energet- 
ically bound to that group; and (4) reject groups that 

^ [http : //«w«-hpcc . astro . Washington . edu/t ools/skid ■ html | 
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have fewer than a minimum number of members. We 
used a 16-member cutoff when running SKID but apphed 
a much more conservative cut at 64 star particles later 
on (§|31). We employ SKID rather than another group 
finder because test runs suggest that it separates satel- 
lite galaxies from main galaxies more accurately than 
halo finders that select objects based on dark matter 
al one. As a check we co mpared our results to those 
of INaeamine et all l|2004ar . who used a different group 
finder with the same simulations, and found no signifi- 
cant differences (§13). 

2.3. Simulated Galaxy SEDs 

For each SKID-identified galaxy, we treat each star 
particle as a single stellar population (SSP) formed at 
the time the star particle was spawned. We then de- 
termine the broadband photometric properties of each 
galaxy by processing its star formation history through 
the GALAXE V library of evolutionary population syn- 
thesis models (|Bruzual fc Ch ariot" 2003V redshifting the 
resulting spectrum, and convolving it with our photo- 
metric filters. We use a'Chabrier ( 2003) IMF. The filters 
we use here are the ACS F435W, F606W, F775W, and 
F850LP (hereafter, ^435, Vgoe, *775: and zgso, respec- 
tively); the VLT /ISAAC J, H, and K,; and the IRAC 
3.6fi, A.bfi, 5.8//, and S.Ofi ([3.6], [4.5], [5.8], and [8.0], 
respectively) bands. Throug hout this work we use AB 
mag nitudes (|Oke fc Gunnlll983iV 

As the GALAXEV library covers a range of metallic- 
ities, we must decide which one to use. Figure ^ shows 
that simulated galaxies' metallicities log(Z/Zo) a,t z — A 
fall between —0.7 and —0.4. Tests showed that using the 
log(Z/ZQ) = —0.7 models for all galaxies produces er- 
rors of less than 0.1 magnitudes for all filters of interest, 
which is good enough to preserve any important quali- 
tative trends in the data. Therefore, we present results 
only from models with log(Z/Z0) = —0.7. Since much 
of the population synthesis work in the literature uses 
models with solar metallicity, we note that using the so- 
lar metallicity models would redden our galaxies by less 
than 0.01 magnitudes in Ks - [3.6] (w rest-frame B — V) 
and dim them by 0.2 and 0.3 magnitudes in rest-frame 
optical (observed [3.6]) and UV (observed ^775), respec- 
tively. 

2.4. Reddening from Dust and the IGM 

The most common method for inferring the star for- 
mation rates of galaxies involves measuring the rest- 
frame UV light emitted by short-lived massive stars. 
Unfortunately, evidence suggests that much of the rest- 
frame UV light from star-forming galaxies is absorbed 
by a dusty interstella r medium and re-ernitted at far- 
infrared wavelengths (fAdelbereer fc Steide]|l2000j) . Both 
the amount of dust in each galaxy and the variation of 
the extinction as a function of wavelength (the "E{B — 
V) distribution" and "reddening law," respectively) are 
poorly constrained for LBGs at present. Available ob- 
servational constraints tend to assume a iCalzetti et al.l 
l|2000l) reddening curve, which is parametrized in terms 
of the color excess E{B — V). Several studies have 
derived E{B — V) d istributions using t his law from 
population synthesis (Shaolev et al. 200 Ij) or UV spec- 
tral slope (Adelberger fc Steidel,2000: Ouchi et alJl2004|) 
techniques. Consensus has centered on a median value 



of E{B — U)~ 0.15 with some scatter for LBGs at 2; ~ 3, 
corresponding to a median factor of 3-5 extinction in the 
UV. For consi s tency with past work, we will employ the 
ICalzetti et al.l ll2000f) redden ing law, although we briefly 
consider the 'Chariot fc Fall (200(/) law as well. 

Given the reddening law, a color excess E{B — V) must 
be chosen for each galaxy. As argued above, it is not ap- 
propriate to simply apply an observed E( B — V) distribu- 
tion such as that of Shaulev et alJ l)200lD to our sample, 
as there is no a priori reason to expect that the intrinsic 
and observed E{B — V) distributions are the same. The 
intrinsic E{B — V) distribution could, for example, in- 
clude a significant population of highly reddened objects; 
the dust in such objects would render their observable 
rest-frame UV colors too red for them to pass the LBG 
color cuts that have been use d to date (see e.g. the BzK 
samples of lDaddi et al]l2004|) . 

The state-of-the-art in intri nsic E(B — V) distribu- 
tions is more art than science. iNagamine et a l. (2004^ 
do not assume an intrinsic E(B — V) distribution, but 
rather plot their results for several different uniform val- 
ues of E{B — V) and note that the data match the 
simulated LFs and number counts well if the typical 
E(B— V) at redshifts z ^ 3-5 is 0.15. By contr ast, recent 
SAMs jSomerville et alJl2001tirdzi et alJ l2004') scale red- 
dening with UV luminosity in a way that mimics the cor- 
relations observed in nearby starburst galaxies; we will 
explore this E{B — V) distribution, but we note that it 
is not well-motivated for our sample because most of our 
simulated LBGs are not starbursts. 

Here, we employ the novel approach of inferring each 
galaxy's dust content E{B — V) from its metallicity. It is 
natural to expect a trend between reddening and metal- 
licity since dust grains are made of metals; indeed, such 
a trend was first observed in lUE spectra of starbursts 
by Storchi-Bergmann et al. (1994^. We use a correla- 
tion observed between the reddening and metallicity in 
the Sloan Digital Sky Survey (SDSS) main galaxy sample 
(C. Tremonti, private communication), extrapolating the 
mean trend to each galaxy's metallicity. For each galaxy, 
we also add a Gaussian scatter SE with variance equal 
to one half of the mean color excess {E{B — V)) at that 
galaxy's metallicity; this produces a good match to scat- 
ter observed in the SDSS sample. The final relation is 



E{B - U) = 9.0 X Z' 



0.9 



SE, 



(1) 



where Z is the galaxy's mean stellar metallicity (ex- 
pressed as a ratio of the mass in metals to the total stellar 
mass) and the scatter SE is of a Gaussian form given by 



fiSE) diSE) = 



1 



UE{B-V))^ 



X exp 



E{B - V) 



^(E{B 



(2) 



Figure n] shows the results of this prescription applied 
to all simulated galaxies by giving the average metallicity 
of the star particles in each simulated galaxy versus the 
applied color excess E{B — V). The solar abundance is 
marked with a dotted line. While the simulation does 
not trace the abundances of individual metal species, we 
report the metallicity as an oxygen abundance since this 
is a common tracer of other metals. 
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Fig. 1. — The applied E (B — V) distr ibution. The x-axis gives the metallicity of each galaxy while the y-axis gives the applied dust 
reddening using the Calzctt i et al.l 120001) reddening curve. Solar metallicity is indicated by the dotted line. Red circles and blue crosses 
are from the G6 and D5 simulations, respectively. Point size scales linearly with log(A/, /Mq) as indicated. Note that more massive objects 
suffer significantly more dust reddening in this prescription. 



As we will show, this method yields results that are 
in broad agreement with other methods, such as ap- 
plying a uniform E{B — V) to all galaxies or picking 
E{B — V) from a Gaussian distribution. It naturally 
produces a weak correlation between mass and redden- 
ing since galaxies with higher stellar mass tend to have 
higher metallicities. In principle, it also has the potential 
to produce massive, evolved galaxies with no gas and (un- 
physically) significant dust. However, inspection of our 
simulated B-dropout sample revealed no massive galaxies 
with abnormally low gas surface density or star forma- 
tion rate. The lack of massive gas-free galaxies at low 
rcdshifts would certainly indicate a failing of the simula- 
tion, but it is unclear whether it is a failing at z ^ 4. In 
any case, it does remove the necessity of further refining 
our reddening prescription at this point. 

Figure 121 shows our fiducial E{B — V) distribution and 



co mpares it with t he me an o f the distr i butio n used 
in iSomerville et alJ l|2001') and ' Idzi et all ^2004) . The 

abscissa gives the observed flux in ^775 before redden- 
ing is applied and the ordinate gives the applied red- 
deni ng E(B — V). T o generate the curve, we adapted 
the ISomerville et al.l (pOOl) prescription assuming a 
;^-o.7 reddening curve and Ry = Av/E{B - V) ^ 
4.05 (Calzetti et al. 2000) to derive E{B -V) = 0.189 x 
-[^00.2(24.98-1775) _ Using the mean relation between the 
unextinguished flux in ^775 and the intrinsic star for- 
mation rate for our galaxies that we derive as de- 
scribed in § 12.31 this can also be written as E{B — V) = 

0.189 X 10(-0-44+0.521og(Af.)) ^j^g^.^ jg ^^-l ^ 

While our prescription naturally leads to UV-brighter 
objects having greater extinction, our mean trend is not 
nearly as strong as theirs. 
For comparison, we explore several other E{B — V) dis- 
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Fig. 2. — Applied reddening E{B — V) versus unextinguished flux 
in *775 (rest-frame 1500 A) for all resolved galaxies. Open triangles 
give our fiducial prescripti on and the solid line g ives the mean 
of the prescription used in 'Somerville et al." f2001'J . Our applied 
reddening scales much less strongly with star formation rate than 
the prescription used in lSomerville et al.i 1.200 L) . 



tributions. In one, each galaxy's reddening is determined 
randomly via a Gaussian distribution oi E{B — V) with 
mean and standard deviation 0.15 but truncated at 0.0 
so that no galaxy suffers negative extinction^; this is re- 
ferred to hereafter as the "random" E{B — V) distribu- 
tion and is similar to the distribution used bv Nigh t et al.l 
((2005). We also use a sample in which each galaxy's red- 
dening is set to 0.12 (the "flat" E{B - V) distribution); 
this number was chosen to produce reasonable agreement 
between the simulated and observed luminosity densities 
in the rest- frame UV (cf. § \A.2\i . Unless otherwise noted, 
we employ our fiducial E{B — V) distribution for all of 
our results. 

Finally, in addition to dust extinction we must model 
attenuation owing to the IGM along the line of sight to 
the galaxy. We use the prescription for the mean at- 
tenuation given by Madau (1995), which accounts for a 
stochastic distribution of optically thin Lyman-a forest 
clouds as well as optically thick Lyman limit systems. 
For our z — 4 galaxies, this prescription suppresses the 
simulated fluxes in the observed -B435 and Veoe bands by 
roughly 1.6 ± 0.1 and 0.4 ± 0.1 magnitudes, respectively 
(the scatter owes to the variation in the shapes of the 
intrinsic galaxy SEDS), and does not affect redder bands 
at all. 

3. NUMERICAL RESOLUTION EFFECTS 

There are a number of ways in which cosmological N- 
body simulations can fail to resolve the physical prop- 
ert ies of galaxies. As a releva nt example, Figure 10 
in iSpringel fc HernauistI (|2£)03b) shows how the simu- 
lated cosmic star formation rate density at a given red- 
shift varies with mass resolution. Broadly, using higher 

^ Note, however, that some authors (e.g., ISteidel et ahl 119991) 
have derived apparently negative values of E{B — V) owing to 
significant Lyman alpha emission boosting the rest-frame UV flux. 



mass resolution allows lower-mass objects to be resolved, 
and allows structure formation at all scales to be resolved 
earlier in a simulation; by contrast, simulating a larger 
cosmological volume typically entails using lower mass 
resolution but allows rare, massive objects to be simu- 
lated. For these reasons, the cosmic star formation rate 
density in the G6 simulation lies below the correct value 
inferred from higher resolution simulations until slightly 
after z = 4. During this time, some of its galaxies may 
possess lower stellar mass and bluer colors than the con- 
verged values. Since we are mainly considering the prop- 
erties of the largest galaxies that collapsed earliest, we 
expect our star formation rates to show better conver- 
gence than the overall galaxy population. However, as 
a check that our simulated sample's properties are re- 
solved, we repeat our analysis with th e higher-resolution 
D5 simulation of ISpringel fc Hernau ist (2003b}, whose 
cosmic star formation rate density converges to the re- 
solved value before z = 4. D5 was started aX z — 159 with 
324^ dark matter and 324^ SPH particles in a cube of co- 
moving side length 33.75/i~^Mpc. It has a SPH particle 
mass resolution of 1.26 x 10^ ft.~^M0 and an equivalent 
Plummer softening length of of 4.17/i~^kpc, and uses the 
same cosmology as G6. 

Mass Resolution: We determined the simulated mass 
resolution by comparing the stellar mass functions for 
the G6 and D5 simulations (Figure |2| top panel) . It 
is tempting to assume that each simulation resolves the 
properties of all galaxies whose stellar mass exceeds the 
mass at which the simulated stellar mass function tracks 
that of the higher resolution simulation; Figure |31 in this 
case tells us that the G6 resolves galaxies with stellar 
mass above log(M:,/M0) = 8.5. Unfortunately, since 
we are interested in star formation histories of galaxies, 
this criterion is insufficient. The bottom panels of Fig- 
ure 13 plot the total baryonic mass versus stellar mass 
for the two simulations. Inspection shows that the cor- 
relation between total baryonic and stellar mass is tight 
at high stellar masses but becomes noisy at lower stel- 
lar masses, well before the stellar mass function turns 
over. This stochasticity, arising as a result of the stochas- 
tic prescription for star formation in the code, would 
result in unphysical scatter in our simulated star for- 
mation histories. To avoid this, we impose a stricter 
minimum stellar mass cut for G6 at 64 star particles or 
log(M*/M0) = 9.64 (the sohd hue in the plot), which 
is slightly below th e stellar mas s obser ved for typical 
(L*) LBGs at z ~ 3 (|Shaplev et al..,.2001:.Papovich et alJ 
12001(1 . Above this limit, the variance in the baryonic- 
stellar mass relation is reasonably small, so it is mostly 
free of stochasticity effects. We test this more thoroughly 
below. 

To fully probe to GOODS depth, we include results 
from the D5 simulation in order to flU out the low end 
of the mass function; using the 64 star particle cri- 
terion on D5 allows us to consider galaxies down to 
log(M*/M0) > 8.76. The normalization of the G6 mass 
function is lower than that of the D5 by about 40% at 
stellar masses below log(M»/M0) = 10.3, likely the re- 
sult of overmerging owing to the lower spatial resolution 
in G6. 

Color resolution: A galaxy's color is determined mainly 
by when its stars formed. In a lower resolution simula- 
tion, a galaxy remains unresolved for longer, so that its 
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Fig. 3. — (top) The mass functions of all galaxies found by SKID in the G6 and D5 simulations, (middle and bottom) Total baryonic 
mass versus stellar mass for all simulated galaxies in the D5 (middle) and G6 (bottom) simulations. The discretized behavior at low stellar 
masses results from finite mass resolution. The heavy vertical lines indicate where our 64 star particle mass resolution cut falls in each 
simulation. Our stellar mass cut removes galaxies whose star formation histories are not well-resolved owing to the stochasticity in the star 
formation prescription. 



star formation is delayed. The galaxy will then be bluer 
than it should be, and may even experience an unphysi- 
cal burst of star formation just as the density threshold 
for star formation is exceeded (M. Fardal, private com- 
munication). Indeed, we find that globally at z = 4, G6 
galaxies have roughly twice the birthrate (i.e. star for- 
mation per unit stellar mass) as D5 galaxies, indicating 
that for the general galaxy population in G6, this effect 
could be significant. 

Does this offset in birthrates propagate into GOODS- 
observable quantities? One measure of the birthrate is 
the Veoe — [3-6] color, which probes the relative amount 
of light from young, OB-type stars to later-type (A and 
later) stars, thus constraining the ratio of the amount 
of recent star formation to the amount from previous 
events. Thus, to answer our question we constructed his- 



tograms of reddened galaxy colors in observed Vgoe ~ [3-6] 
for the two simulations, shown in Figure 01 The plot 
shows only those galaxies whose observed IRAC 3.6/i 
fluxes lie in the range 22 < [3.6] < 24.25 (corre- 
sponding roughly to stellar masses in the range 9.7 < 
log(A4/M0) < 10.5) and which satisfy the GOODS 
z — A LEG color and magnitude cuts (cf. §0)). There may 
be a slight excess at the blue end of the G6 histogram, 
although this may result from an underpopulated tail for 
the D5 distribution owing to its small volume. Overall, 
however, there is remarkable agreement given the dif- 
ferent birthrates; a Kolmogorov-Smirnov test of the two 
color distributions yields a 99.2% probability that they 
are drawn from the same distribution. This suggests that 
the G6 simulated birthrates and colors are generally in- 
sensitive to resolution effects for our choice of a 64-star 
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Fig. 4. — Histograms of Vqo6 — [3.6] colors of simulated B-dropouts with optical magnitudes in the range 22 < [3.6] < 24.25 from G6 (solid, 
red) and D5 (dotted, blue), normalized to unit area. Errors are from jackknife resampling 14.21 . Dust reddening has been applied via 
our fiducial metallicity-derived prescription. The agreement between galaxies' colors in the high-resolution and low-resolution simulations 
indicates that our stellar mass cut removes galaxies whose star formation histories are not resolved in the low-resolution simulation. 



particle cut and GOODS selection. 
Group Finder: The choice of group finder used to iden- 
tify galaxies in the simulated data could in principle af- 
fect the overall properties of the selected sample. Typi- 
cally, however, the baryonic clumps representing galaxies 
are much better separated than dark matter halos in such 
a way that unique identification is not difficult. We com- 
pared LBG number cou nts in the D5 and G6 simulations 
at z = 3 to the results of iNagamine et aTl (2004a b), who 
used an entirely different group finder on the same simu- 
lations. In the G6 and D5 simulations with no reddening, 
they found 12,402 and 202 galaxies that consisted of at 
least 32 particles and satisfied the Steidel et al. (2003) 
color cuts, yielding comoving source densities of -1.91 
and -2.28 in units log(N/ft,^'^]V[pc^). For the same sam- 
ple definition and using our group finder, we extracted 
11,626 and 244 LBGs, corresponding to source densities 



of -1.93 and -2.19. These differences are not easy to in- 
terpret, but as they are small we proceed while bearing 
in mind that the choice of group finder may introduce 
a « 10% uncertainty into volume-averaged predictions; 
this uncertainty is small compared to, for instance, un- 
certainties in our reddening prescription. 

4. PROPERTIES OF THE Z = 4 SAMPLE 

4.1. Sample Definition 

In order to lend concreteness and relevance to our 
simulated sample, we study those galaxies that satisfy 
the GOO DS z ^ 4 LBG color selection criteria given 
in Giavalis coet all l|2004bl) : 

(B435 - Veoe) > 1.2 + 1.4 x (Vgoe - ^sso) A 

(S435 - ^606) > 1-2 A (1^606 - ^850) < 1.2. 

Unless otherwise stated, we combine these color cuts with 
a magnitude cut at zgso < 26.5, which is roughly the lOcr 
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detect ion limit for point sources in 2:350 (jGiavalisco et all 
I2f)n4aj) . Additionally, unless otherwise stated we ap- 
ply dust reddening to our sample using our fiducial 
metallicity-derived E{B — V) distribution. In this sec- 
tion we discuss how these criteria confine the properties 
of the sample under investigation. 

The B-dropout color cuts require that LBGs be red in 
B435 — Vgoe and blue in Vgoe — -^sso- Inspection revealed 
that all of the simulated galaxies satisfy th e redness cut in 
S435 — Veoe- This results directly from the lMadaul l)1995fl 
prescription for attenuation by the IGM, so it is not a 
direct prediction of our simulation. From now on, we 
ignore the B435 — Vgoe color selection criterion. The up- 
per panel of Figure |31 shows the simulated sample, before 
applying color and detection cuts, in the Vgoe — ^sso vs. 
0850 color-magnitude space. The point sizes are scaled 
by stellar mass while the G6 galaxies are red and the 
D5 galaxies are blue. Heavy lines indicate the GOODS 
cuts. The magnitude cut excludes most of the low-mass 
galaxies and a few of the massive ones while the color 
cut independently excludes only a few massive galaxies. 
The slope of the discontinuity between the G6 and D5 
samples indicates the line of constant stellar mass — at 
a given stellar mass, galaxies that are brighter are also 
bluer since both effects result from an increased star for- 
mation rate or a younger stellar population. The bottom 
panel plots the LF for each simulation in units of log(N 
mag~^ Mpc^"^). The smooth connection between the 
two samples where the G6 turns over owing to our stel- 
lar mass minimum gives further evidence that the G6 
galaxies are resolved down to that mass limit. 

Figure shows the mass selection function of the 
simulated GOODS sample for the fiducial and random 
E{B — V) distributions as well as a dust-free case (top, 
middle, and bottom panels, respectively); thus, the top 
panel is a different way of looking at Figure El Since 
none of the simulated galaxies are excluded by the cut 
in — Veoe, we consider only what fraction of the 

galaxies, as a function of stellar mass, are bluer than 
1.2 in Vgoe — -^850 (dotted blue line), brighter than 26.5 
in ^850 (magenta short-dashed line), or both bright and 
blue — that is, observable (solid black line). Not surpris- 
ingly, low-mass galaxies are missed owing to faintness 
while some massive galaxies are missed owing to redness. 
Interestingly, both of these reasonable E{B — V) distri- 
butions predict that the selection functions reach their 
maximum around 10^*^ Mq with a fairly steep dropoff 
to lower masses. For the fiducial case, there is addition- 
ally a shallow decline toward higher masses owing to the 
higher reddening experienced by more massive, metal- 
rich galaxies. In either case, 10-20% of galaxies with 
stellar mass larger than 10^° M© are missed because dust 
causes them to be too red in Vgoe — zs5o- This effect ac- 
counts for much of the difference between reddened and 
unreddened rest-frame UV luminosity functions at the 
bright end of Figure |S1 

The red long-dashed line shows, as a function of stellar 
mass, the cumulative fraction of stellar mass density in 
galaxies at or above a given stellar mass that is selected 
by the GOODS color and magnitude cuts and is in prin- 
ciple observable. This can be thought of as a normalized 
convolution of the mass function 0(to) with the selection 
function /obs(M*): 



P*,obs(> M,) ^ JMr''^'l'Mfobsim)m dm 

Pa.{> M,) jMniax j^j^ 

Several points are of interest. First, the simulations 
predict that, even in an ideal universe with no dust, no 
more than 60% of the total stellar mass density resid- 
ing in galaxies more massive than the resolution limit at 
log(M*/M0) = 8.76 would be observable owing largely 
to low-mass galaxies being too faint. In fact, this is an 
upper limit on the actual completeness over all galaxies 
since the simulated galaxy mass function is likely prolific 
at masses below the resolution limit. Second, once dust 
is added the observed fraction falls to between 30% and 
40%, depending on the E{B — V) distribution. Third, 
some fraction of real galaxies is always unobservable ow- 
ing to blending with stars or other galaxies, further sup- 
pressing the fraction of stellar mass density observable; 
this effect is routinely accounted for via Monte Carlo sim- 
ulations. Finally, it is interesting to repeat the exercise 
for the cumulative star formation density by simply re- 
placing TO with 771 in the above integral, and comparing. 
The total cumulative fractions of stellar mass density ob- 
served are 41, 28, and 60% for the fiducial, random, and 
dust-free reddening distributions, respectively; the total 
cumulative star formation fractions observed are 51, 36, 
and 69% for the same distributions. Thus, the observed 
galaxies may contain a somewhat higher (by « 10%) frac- 
tion of the star formation density than the stellar mass 
density. This is expected since the brightness cut in 0850 
corresponds to a rest-frame UV selection, which is ex- 
pected to trace the star formation rate better than the 
stellar mass. However, the small magnitude of the effect 
indicates that such a selection does not exclude the bulk 
of the stellar mass. 

Previously, iNagamine et al . ( 2004c) compared the to- 
tal density of stellar mass in the G6 simulation to obser- 
vations at z '--^ 3 and found that the simulation includes 
roughly twice as much stellar material as is observed at 
high redshift if the simulation is normalized so as to agree 
with observations at low redshift. They speculated that 
the missed stellar mass may be hidden in massive galax- 
ies that are too red to be observed. Figure El continues 
this line of inquiry. Our findings are in qualitative agree- 
ment in that we show how no more than half of the stellar 
mass is likely to be observable via the Lyman break tech- 
nique. We confirm that a portion of the missed stellar 
mass is hidden in massive, dusty galaxies, although the 
number of such galaxies is sensitive to the amount of 
dust extinction applied. However, as mentioned we also 
find that much — and possibly most — of the missed stel- 
lar mass at this redshift is hidden in numerous low-mass 
galaxies that are simply too faint to be observable rather 
than in massive red or "red and dead" systems. 

It is interesting to compare o ur completeness re sults to 
those of semi-analytic models. Ildzi et all ll2004|) used a 
semi-analytic model to predict that the B-dropout selec- 
tion criteria select less than 50% of the stellar mass den- 
sity in galaxies more massive than log(M*/M0) > 9.57. 
For the same range of stellar mass, we find that the B- 
dropout criteria select 82% and 58% of the stellar mass 
density for the cases of fiducial and random reddening, 
respectively. The difference almost certainly arises from 
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Fig. 5. — (top) GOODS color-magnitude diagram of resolved G6 (red circles) and D5 (blue crosses) galaxies at z=4. Point sizes scale 
linearly with log(M»/M0) as indicated. The heavy lines indicate the B-dropout color- and brightness cuts, (bottom) The zgso hFs for 
the G6 (solid) and D5 (dotted) simulations in units of log(N mag~^ Mpc~^). The B-dropout color and magnitude cuts miss low-mass 
galaxies owing to faintness and massive galaxies owing to redness. 



the very different ways in which the two models treat 
star formation. For example, our LBGs are generally 
more massive than those of Idzi et a l. (2004), and the 
trend between stellar mass and rest-frame UV flux is 
much stronger in our sample than in theirs (cf. their Fig- 
ure 2). Our results qualitatively agree with theirs in the 
sense that most of the missed stellar mass results from 
the magnitude limit at zgso < 26.5 rather than the color 
cuts. 

The fact that a reasonable fraction of massive galaxies 
is missed by the B-dropout color cuts raises the ques- 
tion of whether different selection criteria could isolate a 
larger — or at least different jgortion of the high-redshift 
galaxy population. IDaddi et al.1 l|2004) proposed a two- 
color criterion for isolating high-redshift galaxies that is 
intended to include highly obscured star-forming galaxies 
as well as massive galaxies with old stellar populations. 



Their criterion is RJL = ( J - [3.6]) - 1.4(i? - J) > 0. 
This color cut excludes none of the resolved galaxies in 
our simulations for redshifts z < 4.7 (by z « 4.7, the 
IGM substantially suppresses flux in observed R, ren- 
dering the galaxies too red in observed R — J). Hence, 
our simulations suggest that this may be an effective way 
to select galaxies that are missed by the Lyman dropout 
technique as long as samples are not heavily contami- 
nated by low-redshift interlopers. 

Finally, we remark on the expected sample complete- 
ness in the IRAC bands. While the rest-frame UV data 
probe to zs5o — 26.5, the IRAC data set achieves rough 
5cr detection limits of 26.1, 25.5, 23.5, and 23.4 in the 
[3.6], [4.5], [5.8], and [8.0] channels, respectively. Our 
simulations indicate that with these limits, the fraction 
of B-dropouts that will be detected is 96%, 82%, 12%, 
and 11% for the same bands. Thus, only the [3.6] and 
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Fig. 6. — The mass selection function of LBGs in the simulated GOODS sample for three different E{B — V) distributions: fiducial 
(top), random (middle) and dust-free (bottom). In each plot the blue, dotted line gives the fraction of galaxies in that mass bin that are 
bluer than Veoe — ^850 = 1-2; the magenta, short-dashed line gives the fraction that are brighter than 2850 = 26.5; the black, solid line 
gives the fraction that are both bright and blue; and the red, long-dashed line gives the fraction of stellar mass density in galaxies at or 
above that stellar mass that would be observable in GOODS. For any reasonable intrinsic E{B — V) distribution, the B-dropout selection 
criteria remove low mass (log(Mt/MQ) < 9.5) galaxies owing to faintness and massive galaxies owing to dust reddening. 



[4.5] channels reach deep enough to observe a majority 
of galaxies selected via the B-dropout criteria. 

4.2. Luminosity Density of the Universe 

With the expansive baseline in wavelength afforded by 
GOODS, we can construct a coarse average spectrum of 
LBGs by summing the flux densities of all the galaxies 
through each photometric band. The resulting SED is 
known as the "luminosity density" p^, and can be used 
to study the average properties of th e galaxy population. 
For example. iPapovich et al.l l)2004D used this method to 
suggest that the universe's stellar mass density increases 
by 33 % between z = 4 and z = 3, and Rudnick ct al. 
(|2003'1 used it to infer that the density of stellar mass in 
massive galaxies grows by 10 x between z = 3 and z = 0. 

We compute p^, following the method of lPapovich et all 



lj2004[l . We consider only those galaxies that satisfy the 
B-dropout color cuts, the brightness cut at zsso < 26.5, 
and an additional brightness cut at 1775 < m*(«775) -I- 1, 
wh ere m*(?77s) = 2 5.74 i s derived from the value used 
in IPapovich et al.l i2004l w hich in turn was derived 
from lGia^f!5Co''etalTl20Q4b(l by redshifting from z = 3.9 
to z = 4. The flux densities of the selected galaxies in all 
relevant photometric bands are converted to luminosities, 
added, and divided by the simulation volume to produce 
the luminosity density p^. 

We compared the luminosity densities of the two simu- 
lations computed in this way as a further resolution check 
on the G6 simulation, using only galaxies whose stellar 
masses were in the range 9.7 < log(M*/MQ) < 10.5 since 
this range is resolved by both simulations (cf. Figure |3J). 
We found that, for these galaxies, the D5 luminosity den- 
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Fig. 7. — The luminosity density of the universe for B-dropouts at z ~ 4 in the rest-frame UV to optical bands. The solid, red line gives 
the curve assuming no dust reddening; the magenta, long-dashed line assumes our fiducial reddening prescription; the green, dotted 
line assumes a uniform reddening for all galaxies; the blue, short-dashed line assumes a Gaussian distribution of E{B — V) values (see 
text); the black dot-l ong-dashed curv e assumes the Chariot &; Fall (2000) reddening curve; and the black dot-short-dashed curve gives 
the observations from IFapovich et al.l 12004,) . The fiducial reddening prescription produces broad agreement with the observations while 
showing a small but interesting discrepancy in the observed H-band (see text for discussion). 



sity is 81%-94% that of G6 at this redshift for the vari- 
ous bands. This is smaller than the uncertainty in each 
simulation owing to Poisson noise and cosmic variance, 
which wc derive below, so we ignore this offset. Since 
neither of our simulations completely covers the mass 
range of interest, we compute the total luminosity den- 
sity by summing the independent contributions from G6 
galaxies with log(M*/Mo) > 9.7 and D5 galaxies with 
log(Af,/M0) < 9.7. 

We estimate uncertain ty in p,^ using the jack- 
knife method iLuptonI Il993l: iZehavi et all 120021: 
iWeinberg et al.l 1200411 . In general, the uncertainty 
in estimates of volume-averaged cosmological quantities 
such as the luminosity density receives independent 
contributions from Poisson noise and cosmic variance. 
One simple way to account for both effects is to divide 



up the sampled region into smaller subvolumes of com- 
parable size, estimate in each volume, and take the 
var iance in thes e estim ates as the overall uncertainty. 
As IZehavi et al.l l)2002|) have noted, this method tends 
to overestimate the error owing to cosmic variance on 
scales that are larger than the size of the subsamples. 
The jackknife method utilizes subsamples whose sizes 
are similar to that of the full sample. In this method, 
each subsample is obtained by excluding a small region 
(such as one octant of the simulation) from the full 
sample in turn and estimating p,y in the remaining 
sample. The full uncertainty is then readily obtained 
from the formula: 



2/ ^ 



N 



i=l 



(4) 



High-Redshift Galaxies in Hydro Simulations 



13 



where N indicates the number of subsamples. We obtain 
the uncertainty in the p„ for each simulation by comput- 
ing jackknife errors over octant subsamples in this way 
and, when combining G6 and D5, we add the respective 
errors in quadrature (this is justified since the two simu- 
lations are run with different random number seeds and 
therefore provide independent measurements of the same 
quantity) . 

Figure[7|shows the result. The solid red line is obtained 
in a dust-free universe; the purple long-dashed, green 
dotted, and blue short-dashed lines use the C alzetti et al.l 
(j2003) reddening curve with our fiducial, flat, and 
random E{B — V) distributions, respectivel y; the 
black dot-long-dashed curve uses the Chariot fcFalil 
(|2p00f) reddening prescription; and the black dot-short- 
dash ed curve gives the data of ^Pap ovich et al. (, 200 4) . 
Th e iCalzetti et aLl (2000) curve is widely used, but 
the lCharlot fc Falll ((200a) law is not yet as common de- 
spite its sound physical motivation. In its simplified 
form, it uses a A"'''^ reddening curve and attempts to ac- 
count for the dusty environments around young stars by 
assuming that light from stars that are younger than the 
maximum age of their birth clouds (which we take to be 
10^ years) has an optical depth to dust scattering that is 
greater than the optical depth for the older stars by some 
factor (which we take to be 3). For this curve, we took 
the normalization Ay to be such that the "old" stars for 
each galaxy suffer the same extinction in rest-frame V as 
the entire population in our fiducial metallicity-derived 
case; the young stars are then at three times this optical 
depth. 

It is clear that the fiducial reddening prescription pro- 
duces broad agreement with the data from the rest- frame 
UV into the rest-frame optical. It is also not too differ- 
ent from the results of the other E{B — V) distributions. 
In addition, for reasonable amounts of dust extinction, 
the shape of is not sensitive to the reddening prescrip- 
tion: all of the simulated curves show relatively flat UV 
and optical continua with a strong break at 4000 A. We 
emphasize that the agreement between the data and the 
Pi, from our fiducial E{B — V) distribution has not been 
enforced by hand; the fiducial distribution has only been 
calibrated against nearby late-type galaxies from SDSS. 

While all of the simulated data points are affected by 
dust extinction, the points at 1200 A (observed Vgoe) are 
suppressed by an additional 0.4±0.1 magnitudes (w 25%) 
owing to IGM absorption, as noted in S 12.41 Meanwhile, 
the simulated points at 7200 and 9000 A (observed [3.6] 
and [4.5]) are determined primarily by the stellar masses 
of the simulated LBGs with little dependence on dust 
attenuation. They are reasonably robust, and straight- 
forwardly testable, predictions of our simulation. 

It is interesting that the simulated luminosity den- 
sity curve has a strong break at 4000 A that is not visible 
in the data. The 4000 A break results from a large num- 
ber of metal absorption features at wavelengths A < 4000 
A acting in concert with the Balmer break at 3646 A. As 
it occurs only in the spectra of relatively cool stars, it 
is visible only in stellar populations older than ^ 100 
Myr and is widely used as a characteristic age indicator. 
The stellar populations in our galaxies have character- 
istic ages of 300 Myr (Figure lll|l ; therefore a break is 
expected. Although star formation histories of LBGs are 
currently poorly constrained by observations, the work 



of lShaplev et alJ l)200H) suggests that a substantial popu- 
lation of galaxies whose mean stellar ages are old enough 
to show a 4000 A break should be visible by z ~ 4; thus, 
one might expect to see this feature in the data. The in- 
gredients that are most likely to give rise to the discrep- 
ancy are redshift scatter, dust reddening, the simulated 
star formation histories, and uncertainty in the data. As 
the 4000 A break is an important feature for character- 
izing stellar populations, we consider these possibilities 
in turn. 

Redshift scatter is present in the data since the red- 
shifts of LBGs are known only to within the width of 
the selection function resulting from the color cuts that 
are used to define the LBG sample. We find that the 
B-dropout color cuts select most galaxies with 3.5 < z < 
4.5 and very few galaxies outside this range; thus, the op- 
tical data in Figure can be viewed roughly as a convo- 
lution of the unknown, correct luminosity density curve 
with a tophat function of width 400-900 A, depending 
on the observed band. We tested whether this could re- 
move the 4000 A break by introducing artificial scatter 
of ±0.5 into the simulated galaxies' redshifts before ap- 
plying the B-dropout selection criteria and found that 
the simulated luminosity density dropped in normaliza- 
tion but did not change shape. The normalization drop 
results from moving galaxies into redshift ranges where 
the selection function is lower, and the 4000 A break is 
unchanged because the width of the selection function is 
not broad enough to allow significant flux from redwards 
of the break to spill into the observed H band (rest- frame 
3200 A). 

As always, dust reddening is a possibility. With enough 
dust any galaxy's observed UV continuum can be red- 
dened to the point that the 4000 A break disappears. 
Unfortunate ly, if we add enough dust to our galaxies 
(using the ICalzetti et aL .2000 law with a reasonable 
intrinsic E{B — V) distribution) that the shape of the 
simulated luminosity density curve approaches the ob- 
served shape, its normalization drops more than an order 
of magnitude below the observed normalization owing to 
too many galaxies being reddened out of the observed 
sample. Thus for this explanation to be viable, the real 
LBGs must be significantly more luminous and dustier 
than our simulated sample; observations in the [3.6] and 
[4.5] bands should constrain this possibility. An alterna- 
tive check on the normalization is provided by the rest- 
frame UV luminosity function fS I4.3(l . which shows broad 
agreement given our fiducial E(B — V) distribution. Fig- 
ure 13 thus tells us that the simulated 4000 A break is 
insensitive to reddening over a variety of E{B — V) dis- 
tributions and two well-calibrated reddening curves if we 
wish to preserve agreement between the observed and 
simulated rest-frame UV luminosity functions. Hence, 
we do not think that reddening is the main cause al- 
though it may contribute. 

Our currently favored explanation for the "miss- 
in g" 4000 A break is the uncertainty in the data 
of iPapovich et all l|2004D . Their sample was culled from 

50 square arcminutes of data (or « 5% of the comov- 
ing volume in our simulation) owing to the limited areal 
coverage of the VLT /ISAAC data, and their bootstrap 
uncertainties do not account for cosmic variance (note, 
however, that the break was not observed in their U- 



14 



Finlator et al. 



dropouts at 2: '--^ 3 either) . More insight into this problem 
could be gained if the luminosity density measurement 
were reproduced with the full GOODS data set includ- 
ing deep imaging in the critical VLT /ISAAC H-band. 

4.3. Luminosity Functions 

In this section we present rest-frame UV and opti- 
cal luminosity functions for the simulated B-dropout 
sample. Several studies of the luminosity functions of 
LBGs in these simulations have already been published. 
iNagaminc ct al. ( 2004a Ij) found that the simulated lu- 
minosity functions at z ~ 3 in observed V and R repro- 
duce the observed luminosity function at the bright end 
reasonably well while exhibit ing; st eep slopes a '--^ — 2 at 
the faint end. iNight et al.l l)200 5'l performed a similar 
study for LBGs at z = 4 - 6 and found good agreement 
between simulated and observed lumin osity functions in 
i' and z' provided that iCalzetti et al.l l|200Q) reddening 
was assumed with E{B - V) = 0.15 - 0.30. Addition- 
ally, they found a steep slope a ~ —2 at the faint end of 
the simulated luminosity function. At the risk of some 
overlap with their work, we will present a complete dis- 
cussion of the LBGs in these simulations under a single 
set of assumptions such as our mass cuts and assumed 
E{B — V) distribution. In order to create a single lu- 
minosity function, we combine the results of the D5 and 
G6 simulations following the same method as described 
in S lO 

Reddening: Figure [S] shows the simulated 1775 luminos- 
ity function for a number different E{B — V) distribu- 
tions. The galaxies used in creating this plot satisfy 
the GOODS color cuts but do not necessarily satisfy the 
GOODS detection limit zg5o < 26.5. The solid red line 
shows the luminosity function derived if dust is absent 
so that E{B — V)=Q. The magenta long-dashed, green 
dotted, and blue short-dashed lines show the luminos- 
ity function under the assumptions of a fiducial, flat, 
and random E{B — V) distribution, respectively (§13; 
and the cyan dot-short-dashed line assumes that the red- 
dening varies with star formation rate in the way de- 
scribed by Som ervill e et al. (2001). All of these curves 
use the iCalzetti et al.l (2000) reddening curve once the 
reddening E{B — V) for each galaxy is chosen. The black 
long-dash-dotted line uses the Chariot & Fall ( 2000) red- 
dening curve as described in § 14.21 

We briefly consider the various reddened curves in 
turn. Inspection reveals slight differences between the 
results from the flat, random, and fiducial E{B — V) 
distributions. Our fiducial E{B — V) distribution natu- 
rally produces a weak positive correlation between stellar 
mass and reddening that is absent in the random dis- 
tribution, preferentially reddening massive galaxies out 
of the sample (cf. Figure 01. For the same reason, at 
the faint end the random E{B — V) distribution causes 
more galaxies to be reddened out of the sample than 
the fiducial distribution. The random distribution al- 
lows a few more bright galaxies to have low dust con- 
tent than either the flat or fiducial E{B — V) distri- 
butions. The curve from the flat distribution is simi- 
lar to the others except at the bright end, where, once 
again, the redder intrinsic colors cause these galaxies to 
be preferentially reddened out of the sample if high red- 
dening is enforced. Despite these minor differences, it 
is clear that the random, fiducial, and flat E{B — V) 



distributions produce qualitatively sim i lar re sults. The 
curve derived from the lSomerviUe et all polll E{B-V) 
distribution reddens rapidly star-forming galaxies out of 
our sample far too effectively; it is clearly not appropri- 
ate foi_our simulated galaxies. The curve derived from 
the lCharlot fc Fall (JOM) reddening prescription has al- 
most exactly the same shape as our fiducial curve but 
is fainter by roughly 0.5 magnitude; this is the result 
of adding extra extinction to the youngest stars in each 
galaxy. It would no t be difficult to modify the parame- 
ters we used for the lCharlot &: FalL ( ,200 0) law to obtain 
better agreement, but as little insight would be gained 
from this exercise, we do not do so. 

As with the observed mass fraction, the dominant 
source of uncertainty in these predictions is clearly the 
unknown form of the intrinsic E{B ^V) distribution and 
reddening curve in high-redshift galaxies. However, the 
important point is that the overall shape of the simulated 
luminosity function is fairly insensitive to the details of 
the reddening procedure such that reasonable distribu- 
tions agree to within a factor of 2. This gives us further 
confidence that our fiducial E{B — V) distribution will 
yield useful predictions. 

Comparison to Observations: Figure El shows a compar- 
ison of the simulated 1775 luminosity function using the 
fiducial E{B — V) distribution with various observational 
results. The red crosses with jackknife errors indicate the 
simulated GOODS sample, while the arrow at the faint 
end points to the value obtained by relaxing the bright- 
ness cut zgso < 26.5. The plot also includes observed 
LFs from the FORS survey ((Gabasch et al..,200^ at rest- 
frame 1500 A, the Subaru Deep Field l|Ouchi et alJl200l 
at rest-frame 1700 A, and Steidel et al. (jSteidel et alJ 
1999) in observed /; these are indicated by solid, short- 
dashed, and long-dashed curves, respectively. The FORS 
LF was constructed at z = 4 using the fitting function 
gi ven by Equation 1 a nd the values given by Table 4 
inlGaba sch et all l)2004j) . The Steidel et al. z - 4 LF was 
constructed from the updat ed Schechter pa rameters for 
U-dropouts (z ~ 3) given in lAdelberger fc Stcidcl (200(|) 
and following the method described in Ste idel et all 
(ES99): we used a distance modulus to shift to* from 
z = 3.04 to z = 4.0, multiplied the normalization (jj^ 
by 80%, and left a constant, yielding to*, a) — 
(3.52 X 10-3, 25.02, -1.57). The Subaru D eep Field LF 
is the LF for BRi-LBGs from Table 4 of lOiichi et alJ 
l)2004j) . where we have used their derived faint-end slope 
a = —2.2. Note that because the effective wavelengths 
of the bands used in these surveys as well as the mean 
redshifts of their LBG samples differ, offsets of 0.1-0.3 
magnitudes between their reported LFs do not necessar- 
ily indicate disagreements. 

The observed and theoretical GOODS LFs both flatten 
out for magnitudes fainter than 26. By contrast, the 
magnitude-unlimited GOODS luminosity function seems 
to turn over at a fainter magnitude. In fact, this turnover 
owes to the conservative cut in stellar mass: The more we 
relax this restriction, the fainter the magnitude at which 
the rest-frame UV luminosity function turns over. Hence, 
our simulation is predicting an extremely steep intrinsic 
faint-end slope {a ^ —2) that cannot be observed in B- 
dropout samples of GOODS depth. 

LBGs brighter than 25th in 1775 also seem to be over- 
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Fig. 8. — The rest-frame UV LF of B-dropouts assuming various different E{B — V) distributions. The solid red curve assumes 
no reddening; the dotted green curve assumes a uniform reddening; the short-dashed blue curve assumes a Gaussian distribution of 
E{B — V) values (see text); the long-dashed magenta curve assumes our fiducial reddening prescription; the cyan dot-short-dashed curve 
assumes the E{B — V) distribution of^ Somerville ct al. 1 2001); the black dot-long-dashed curve assumes the Chariot & Fall 1 2000) reddening 
curve. For a variety of reasonable E{B — V) distributions, the simulated LF varies by < 2x from the curve that assumes our fiducial 
E(B - V) distribution. 



produced in the simulation. This is analogous to the 
bright-end excess noted in the same simulation at z = 2 
by Nagaminc ct al. ( 2005a) . The excess could result from 
any of four possibilities: (1) incorrect prescription for 
dust reddening; (2) lack of truncation of star formation, 
owing to e.g. feedback from black hole growth; (3) over- 
merging; and (4) incorrect feedback model for star for- 
mation. We consider each of these possibilities in turn, as 
they provide interesting insights into high-redshift galax- 
ies. 

The first possib i lity is consistent with the idea, noted 
by iShaplev et al.1 l)200lD and others, that dust extinc- 
tion correlates with luminosity, and hence in our model, 
star formation rate. Clearly, if our assumed extinction 
scaled more strongly with star f ormation rate, closer to 
that in iSomerville et al.| l)2001l) . then the discrepancy 



would be alleviated. E vidence in support of this pos- 
sibility is provided by iDaddi et al.l |20^): For their 
actively star-forming BzK sample at z >~ 1.4 with 
K < 20, they find a median star formation rate of 
200 Mq yr~^ and a median reddening of E{B — V) ^ 
0.4. Indeed, iNagamine et alJ l)2005b(i found that the 
simulations can reproduce the comoving number den- 
sity of Extremely Red Objects (EROs) at z = 1- 
2 if a uniform reddening of E{B ~ V) ~ 0.4 is as- 
sumed for the entire simulated sample. Using our 
fiducial E{B — distribution, simulated galaxies at 
z = 4 with star formation rates between 150 and 250 
Mq yr^^ suffer a significantly weaker median redden- 
ing of roughly 0.15. Additionally, we note that some 
of the most massive, rapidly star-forming galaxies could 
correspond to extremely dusty starbursts that are ob- 
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Fig. 9. — (large) The simulated rest-frame UV LF of B-dropouts (assuming our fiducial reddening prescription) compared with several 
published curves and the GOODS data. Red crosses indicate the simulated LF in observed 4775; the arrow points to the value of the 
LF at that point if the brightness cut in ^850 is relaxed. Errors in the simulated plot arc from jackknife resampling. LFs from the 
FORS survey (Gabasch ct al. 2004) at rest-frame 1500 A, the Subaru Deep Field lOuchi ct al. 2004) at rest-frame 1700 A, and Steidel 
et al. tSteidel et a l. 1999) in observed / are indicated by solid, short-dashed, and long-dashcd curves, respectively, and were derived from 
the literature as described in the text. The GOODS B-dropout sample of Giavalisco ct al. 1 2004b) in observed 4775 is indicated by blue 
triangles, (inset) The LFs of the D5 and G6 simulations, which were combined to produce the full curve as described in the text. The 
simulated LF assuming our fiducial E{B — V) distribution produces broad agreement with available constraints, with a possible excess at 
the bright end (see text for discussion). 



served as sub-millimeter sources (iChapman et al.ll200^ . 
as we will consider further in § 15.41 Recent obser- 
vations suggest that the effects of dust on these and 
other massive galaxies (log(M,/M0) > 11) may not 
even follow sim ple foreground screen attenuation laws 
such as that of iCalzetti et al.1 l|200nD (jChanman et al.1 
120051 Papovich et al. 200 5, in preparation). In this case, 
the lCharlot fc Falll 1)2000(1 reddening curve is almost cer- 
tainly more appropriate; unfortunately, observations to 
date have not constrained the parameters needed to use 
it at high redshift. 

The second possibility is that our simulation may 
fail to account for some mechanism that suppresses or 



truncates star formation in real galaxies. In this case, 
the bright-end excess is an early sign of the difficulties 
that are commonly seen in cosmological simulations and 
SAMs (Somcrvillc 2004; Naeamine ct al. 2005a). Up- 
comi ng work incorp orating AGN feedback may help with 
this lini Matteo et a,1. 2005: Springel et al.„..2Q05a.b) . 

The third possibility is overmerging owing to poor spa- 
tial resolution. Figure |3| shows that galaxies with stel- 
lar masses log(M*/MQ) < 10.25 are slightly underrepre- 
sentcd in the G6 simulation relative to the D5. Owing 
to its poorer spatial resolution, galaxies that form in the 
G6 are more extended than in D5, rendering them more 
susceptible to tidal disruption and subsequent accretion 
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by more massive galaxies. Massive galaxies in G6 then 
grow more rapidly than they would in a higher-resolution 
simulation, leading in principle to a brighter character- 
istic magnitude M*. The impact of overmerging on the 
bright end of a rest-frame UV LF can be characterized 
by inspecting the galaxy mass function at the massive 
end. Figure |21 shows that, up to the maximum galaxy 
mass produced in the smaller, higher-resolution D5 sim- 
ulation, the mass functions of the D5 and G6 simula- 
tions agree quite well. Together with the evidence that 
the birthrates in the G6 are numerically resolved (Fig- 
ure El), this gives us confidence that overmerging is not 
the dominant source of the discrepancy. 

The final possibility is an incorrect model for supernova 
feedback. In particular, our simulation's prescription for 
kinetic feedback endows all galactic outflows with a com- 
mon velocity, 484 km s~^. This ma,y be an oversimplifi- 
cation. For example, iMartinI l)2005(i found that outflow 
velocities in a sample of 18 ultraluminous infrared galax- 
ies correlate well with star formation rate, and suggested 
that this may result in winds more effectively suppressing 
star formation in massive galaxies. 

In summary, our study of the rest-frame UV LF con- 
flrms the previously cited studies in the following re- 
spects: (1) the simulated and observed LFs show broad 
agreement when a reasonable amount of dust extinction 
is assumed; (2) the intrinsic faint-end slope of the simu- 
lated LF is ~ —2; and (3) an excess at the bright end of 
the simulated LF suggests that more work is required 
to understand the effects of star formation and AGN 
feedback as well as dust reddening on massive galaxies. 
The fact that our work reproduces these results despite 
our having employed a different group finder and a more 
physically-motivated prescription for dust reddening sug- 
gests that these results are robust to the analysis proce- 
dure. 

Figure irni shows the rest-frame optical luminosity func- 
tions and compares them with the rest-frame UV lumi- 
nosity function. The observed [3.6] and [4.5] luminosity 
functions are very similar, suggesting that measurements 
of stellar mass may readily be made in the more sensi- 
tive [3.6] channel. Since LBGs tend to be 1-2 magni- 
tudes brighter in rest-frame optical than in rest-frame 
UV, the rest-frame optical luminosity function is essen- 
tially the same as the rest-frame UV luminosity function, 
shifted 1-2 magnitudes brighter. This is consistent with 
the trend, discussed below, that LBGs are simply the 
most massive of the early galaxies, and that their star 
formation rates correlate with their masses. A closer look 
suggests that the rest-frame UV luminosity function ap- 
pears to be slightly steeper than its optical counterpart, 
consistent with the tendency for more massive galaxies 
to be redder. The simulations predict that the faint end 
of the intrinsic luminosity function is essentially a power- 
law down to much fainter masses than GOODS can ob- 
serve; the turnover in the intrinsic luminosity function at 
[3.6] « 26.5 owes to the 64 star particle mass resolution 
cut. Similarly, the turnover in the observed luminosity 
function at [3.6] « 25 results directly from the GOODS 
magnitude cut at zs5o < 26.5 and translates into a mass 
cut at log(M*/M0) w 9.5 fFigurelHl see also FigureEJ. 

5. COMPARING PHYSICAL AND PHOTOMETRIC 
PROPERTIES OF LBGS 



In this section we focus on correlations between the 
properties of the simulated sample, comparing to the 
data wherever possible. Unless otherwise noted, we em- 
ploy our fidu cial, metallicity-deriv ed E{B — V) distribu- 
tion with the'Calzetti et al." (2000) reddening curve. We 
occasionally compare the physical properties of our sim- 
ulated sample at z = 4 with the properties derived from 
observations at lower redshift z ~ 3 since LBGs at this 
epoch are well-studied. This may not be a valid compar- 
ison since some evolution may occur in the population 
of star-form ing galaxies between these two epochs (e.g., 
lPaDOvich~et al. 20Q3)- On the whole, how ever, the two 
j)opulations are expected to be very similar l|Steidel et alJ 
1999); thus, we proceed while noting that factor of 2 dif- 
ferences do not necessarily suggest a failing of the simu- 
lation. 

5.1. Physical Properties 

Figure 1111 shows a number of physical correlations in 
addition to the basic properties of the simulated GOODS 
z ^ 4 sample. Red points and blue crosses denote G6 and 
D5 galaxies, respectively. The break in the G6 sample 
at log(Af*/M0) = 9.64 owes to our 64 star particle cut. 
The histograms are constructed by considering only G6 
galaxies more massive and D5 galaxies less massive than 
log(M*/M0) = 9.7 and weighting the contribution from 
the D5 galaxies by the ratio of the simulation volumes. 
The contours were constructed using the same method 
and define regions enclosing 63% and 95% of the galaxies. 
For each galaxy, log(M*/MQ) is the sum of the masses of 
the star particles; the star formation rate is the instanta- 
neous star formation rate; is the median age of its star 
particles; 12 -I- log{0/H) is the mean metallicity of its 
star particles; and p/p is the local space density of stel- 
lar mass in B-dropouts (roughly equivalent to an optical 
luminosity- weighted number density) , normalized by the 
mean over the simulation volume. We compared instan- 
taneous star formation rates with star formation rates 
obtained by averaging over the last 100 Myr and found 
that they arc within a factor of 2 15.4(1 . in agreement 
with Weinberg et al. ( 2002). 

Stellar Mass: The stellar masses of our sample fall pre- 
dominantly in the range log(M<,/M0) = 9-10 with a 
peak in the range 9.5-9.8 and a significant tail out to 
log(M,/M0) « 11. The falloff at lower masses owes en- 
tirely to the brightness cut at zs5a < 26.5 (cf. Figures [3 
andO. Overall, this range is slightly more m assive than 
the sa mple from the semi-analytic model of lldzi et all 
(|2004 Figure 2); in particular, those se mi-analytic mod- 
els do not reproduce our high-mass tail. iPapovich et alJ 
pool) used population synthesis techniques to infer the 
properties of the stellar populations of 33 LBGs in the 
HDF-N between z 2 and z ^ 3.5 and found a typi- 
cal steflar mass of 10'^°Mq. iShaplev et aD l|200H) used 
similar techniques to infer the properties of 74 fairly 
bright LBGs at z ~ 3 and found masses in the range 
log(M*/M0) = 9-11 with a median at 1.2 x IO^^Ii-'^Mq. 
While these measurements are somewhat more massive 
than the median in Figure ^2 we note that both of these 
samples are biased towards brighter and thus more mas- 
sive galax ies than our sim ulated GOODS sample. Also, 
given that IPapovich et all 1(2004^) report a '-^ 33% buildup 
of stellar mass between z = 4 and z = 3, it is reasonable 
to suppose that galaxies at z = 3 might be more some- 
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Fig. 10. — (large) The simulated rest-frame UV and optical LFs of B-dropouts. The arrows point to the values of the LF if the brightness 
cut in 2850 is relaxed. Errors are from jackknife resampling, (inset) The luminosity functions of the D5 and G6 simulations in observed 
[3.6], which were combined to produce the top curve as described in the text. The rest-frame optical LF is very similar to the rest-frame 
optical LF, shifted roughly 1—2 magnitudes brighter. 



what more massive than galaxies at z ~ 4. Thus, we 
conclude that our LBG model seems to be able to match 
the distribution of observed steUar masses. 
Star Formation Rates: The bulk of the observed galax- 
ies are forming stars at rates of 1-10 Mq yr~^. A small 
active population is forming stars at aggressive rates of 
100 ^ 1000 Mq yr~^ with the median instantaneous 
star formation rate of the ten most rapidly star-forming 
galaxies being ~ 450 Mq yr~^. 

Our models produce a tight correlation between stel- 
lar mass and star formation rate, a s noted elsewhere 
in hydrodynamic simulations (e.g.. [ Dave et aP 119991: 
IB'cinbcrg eLal. 2002; Nagamine et al J l2004a^~ A fit to 
this relation for our simulations gives 

log(M,/Mo) « 1.141og(M,/Mo) - 10.2, (5) 

which is shown in the plot. As mentioned in § |31 a 
comparison of the mass functions (Figure I^J and the 



birthrates (not shown) of the D5 and G6 galaxies sug- 
gests that the G6 galaxies may be suffering slightly 
from overmerging, but we believe that the star forma- 
tion rates in the G6 simulation cannot be overpredicted 
by more than a factor of 2, which is comparable to 
the intrinsic scatter. Figure thus predicts that a 
few galaxies with instantaneous star formation rates of 
500-1000 Mq yr-i may faU within the GOODS field. 
According to our prescription for dust reddening, they 
should all be observable as LBGs; however, considering 
observational evi dence that the most massive galaxies 
may not o bey the[C alzetti e t al.l l)200CI|) reddening curve 
ifChaoman et al J 120 05: Papovich 2005, in preparation), 
this may not be the case in reality. We will return to 
these massive star formers in § 15.41 For now, the strong 
correlation between stellar mass and star formation rate 
implies that LBGs observed at z = 4 are predominantly 
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Fig. 11. — Physical properties of B-dropouts versus stellar mass, (top) instantaneous star formation rate; (second) median age of the star 
particles; (third) mean metallicity of the star particles; (fourth) local density of observable stellar mass, relative to the mean. Red points 
and blue crosses denote G6 and D5 galaxies, respectively. The contours enclose 63% and 95% of the simulated galaxies. The histograms 
and contours are constructed by volume-weighting the contributions of the D5 and G6 samples as described in the text. Note the strong 
correlation between star formation rate and stellar mass. 



massive galaxies. 

The scatter in the plot of instantaneous star formation 
rate versus stellar mass can be regarded as a measure 
of the star formation rate "burstiness" as a function of 
stellar mass. For example, galaxies with stellar mass of 
10^° M0 have star formation rates between 10 and 30 
Mq yr^^. This means that such a galaxy will, on av- 
erage, tend to form 10 Mq yr~^ although it will some- 
times experience bursts of up to 30 Mq yr^^; if the sim- 
ulation produced more dramatic bursts than this then 
there would be more spread in the plot. Note, however, 
that the simulation does not adequately resolve merger- 
driven starbursts (e.g. Hernquist & Mihos 1995; Mihos 
& Hernquist 1996; Springel et al. 2005b). For those 
galaxies undergoing mergers, the instantaneous star for- 
mation rate reported by the simulation may be an un- 
derestimate even though the total stellar mass formed 



by the end of the me rger event is reliable. Relatedly, 
iWeinberEf et al.l l|2002f) notes that there could be merger 
events between galaxies not massive enough to be re- 
solved in the simulation, which cause them briefly to 
be luminous enough to satisfy the B-dropout magnitude 
cuts. However, the broad agreement between the simu- 
lated and observed rest-frame UV luminosity functions 
suggests that the bulk of the LEG population can be 
accounted for without resorting to such phenomena. 

The lower end of the distribution of star formation 
rates is reasonably consistent with semi-a nalytic re- 
sults (jSomerville et alJl200Hlldzi et alJl2004j) . However, 
the space density of more rapidly star-forming galaxies 
10-100 Mq yr-i is higher in our results, no doubt owing 
to the drastically different star formation and redden- 
ing prescriptions involved. Observationally inferred LEG 
star formation rates are in reasonable agreement with 
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our results ijShaplev et al.ll200H iPapovich et"al]l200lD 
The h istogram of star formation rates in lshapfev et al.l 
()200H) seems to show a higher space density of galaxies 
with star formation rates greater than 10 M0 yr~^, but 
since this sample was hand-selected to span the range 
of possible properties of LBG stellar populations rather 
than to satisfy an unbiased selection criterion, such an 
excess of galaxies exhibiting more extreme properties is 
unsurprising. 

Metallicity: The simulated galaxies' stellar metallic- 
ities correlate tightly with stellar mass. The simu- 
lation computes metallicities by allowing the metal- 
licity of a star-forming gas particle to increase dur- 
ing each timestep in proportion to its current star 
formation rate and a constant yield that is set to 
the solar metallicity 0.02; as mentioned, a star par- 
ticle inherits the metallicity of its parent SPH parti- 
cle. Given that a galactic mass-metallicity relation is 
not assumed in this prescription, it is of interest to 
examine the trend that nonetheless results. Such a 
trend is requir ed by a growing bo d y of observational 
evide nce Ce.g.. iLeaueux et al.l Il97^ lOarnett fc Shieldsl 
119871 iTremonti et all 20041): for a review of this subject 
see the discussions in ITremonti et al. ( 2004). 

The dashed line gives the median relation be- 
tw een stellar mass and g as-phase metallicity observed 
by ITremonti et"al] l)2004|) in their sample of ~ 53, 000 
nearby (z « 0.1) star-forming galaxies. Despite the sim- 
ple prescription for enrichment, the slope of the trend in 
the simulated sample shows excellent agreement with the 
observed low-redshift trend. Over a range of two orders 
of magnitude in stellar mass, the simulated galaxies show 
a stellar metallicity that is roughly 0.6 dex below the low- 
redshift value for their respective stellar masses and the 
most massive objects (log(M*/Mo) > 11) already ex- 
hibit solar metallicity. Closer examination suggests that 
th e simulation may not p roduce the flattening observed 
by ITremonti et a'lTl)2004|) at stellar masses greater than 
log(M*/M0) = 10.5. While such a disagreement may be 
a failing if observed in gas-phase metallicities at low red- 
shift, it is not necessarily a failing for stellar metallicities 
at high redshift. 

Physical processes that could contribute to the ob- 
served mass-metallicity relation include selective re- 
moval of metals by supe rnova-dr iven winds ("deple- 
tio n," first modeled by iL arsonl [1 974 and favored 
by ITremonti et an i2004'l. ISM e nrichment owing to su- 
pernova feedback ("astration," .Leaueux et alJlT97^ . di- 
lution owing to infall of unenriched gas, and late-time re- 
turn of relatively unenriched gas to the ISM from evolved, 
low- mass stars ("recycling"). Of these processes, deple- 
tion contributes in the simulation via the superwind pre- 
scription, astration is built into the enrichment prescrip- 
tion, dilution can contribute, and recycling is not ac- 
counted for. If astration dominated the trend then the 
stellar metallicities in all simulated galaxies would be ex- 
pected to be roughly solar owing to the constant yield 
assumed by the star-formation prescription. If dilution 
dominated the trend then massive galaxies, which are ac- 
creting gas much more rapidly than less massive galaxies, 
would be expected to exhibit lower metallicities — the op- 
posite of what is seen. Thus, the simulation's superwind 
prescription is the most likely cause of the stellar mass- 



metallicity relation, and the excellent agreement between 
the slopes of the observed and simulated relations give us 
further confidence in the superwind prescription as well 
as in the predicted LBG metallicities. 
Median Age: Figure ^2 (second panel) shows that the 
median age of a galaxy's stars is at best weakly correlated 
with stellar m ass. This means th at the phenomenon of 
"downsizing" ijCowie et al. 11199^ . i.e. that larger galax- 
ies are observed to have an older stellar population out 
to z ^ 1, is at best weakly evident in our models at 
z = 4. This is a fundamental and testable prediction of 
our simulations. 

The age distribution extends from 200-600 Myr with a 
median around 350-400 Myr. This is roughly a factor of 
2 older than the populat ion in the semi-analytic simula- 
tions of lldzi et alJ l)2004j) . consistent with the picture in 
which LBGs are massive galaxies rather than starbursts. 
These values are also roughly a factor of two older than 
th e constraints d e rived from observed LBGs at z 3 
by iShaplev et a l.' (2001). However, the uncertainties in 
star formati on histories derived via t his technique are 
substantial l|ShaDlev et al.ll20011 l2?ffll : thus, agreement 
within a factor of two gives us confidence that our models 
have reasonably realistic star formation histories. 
Environment: The fourth panel gives the galaxy den- 
sity p/p in each galaxy's local environment. We define 
the density for each galaxy as the local density (relative 
to the mean) of stellar mass in the 16 nearest galaxies 
that would be observable by GOODS, smoothed with an 
SPH kernel; this is broadly equivalent to a number den- 
sity weighted by rest-frame optical luminosity. There is 
a weak correlation between stellar mass and local matter 
density in the sense that the most massive galaxies are 
in denser environments. Galaxies in denser areas also 
show redder colors and higher star formation rates (Fig- 
ure^)), likely arising as a secondary correlation owing to 
the tight relationship between stellar mass and star for- 
mation rate. Although there is considerable scatter, such 
a correlation might be observable in the GOODS sample. 
We note that weighting by stellar mass (or, equivalently, 
rest-frame optical flux) is critical; when we compute den- 
sities using unweighted number counts these weak trends 
disappear altogether. 

5.2. Physical Versus Photometric Properties 

Figure El compares some physical and photomet- 
ric properties of the simulated sample as it would be 
observed with reddening (solid contours) and without 
(dashed contours). Physical properties include stellar 
mass, star formation rate, age, and local galaxy density. 
Photometric properties include rest-frame UV (Vgoe) 
and optical ([3.6]) magnitudes, a short-baseline rest-UV 
color (Vgoe — ^85o) and a long-baseline UV-optical color 
(^606 — [3-6])- While the GOODS measurements cover 
all four IRAC channels, we plot only the [3.6] band here 
because the [3.6] band will detect the largest fraction of 
the B-dropouts (§01). We note that the trends observed 
in [3.6] persist in the redder bands. Figure IT^ shows the 
following interesting properties: 

- Rest-frame optical data are a better predictor of star 
formation rate (SFR) than rest-frame UV data. This 
seems surprising at first, but is expected given the tight 
stellar mass-star formation rate correlation in Figure ITTj 
along with the fact that reddening introduces much more 
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Fig. 12. — Physical versus photometric properties of B-dropouts. Solid contours enclose 63%, 86%, and 95% of the galaxies assuming 
our fiducial E{B — V) distribution; dotted contours assume a dust-free universe. Contours are constructed by volume-weighting the 
contributions of the D5 and G6 simulations as described in the text. The dashed lines in the plot of stellar mass versus [3.6] flux show 
constant mass-to-light ratio corresponding to 1/8 (brightest), 1/4 (middle), and 1/2 (dimmest) of the solar value. Note that star formation 
rate correlates more strongly with optical flux than with UV flux, even when reddening has not been applied. 



scatter into the UV than into optical light. One may 
wonder whether rest-frame UV would be a better pre- 
dictor of SFR if it were possible to remove the effects of 
reddening accurately; the dotted contours in Figure 1121 
suggest that this would at best allow UV flux to be an 
equivalently good predictor, but not better. 
- The tight correlation between star formation and stel- 
lar mass also gives rise to a correlation between rest- 
frame UV flux and stellar mass, in qualitative agree- 
me nt with the tre n d not e d in U-dropouts at z ~2-3.5 
bv lPapovich et all (|200lD . iShaolev et all l|2nn5j) do not 
find such a correlation in their UV-selected sample at 
z ~ 2. However, Papovich et al. (2005, in prepara- 
tion) find a significant correlation between stellar mass 
and star formation rate in their sample of near-infrared- 
selected galaxies a.t z ^ 1-3.5. Thus, a careful assessment 



of selection biases is probably ne cessary for determining 
whether the lShaplev et al.l 1)2005(1 result is in conflict with 
the simulations. 

- More massive galaxies are redder in both the unred- 
dened and reddened cases, with the trend accentuated 
in the reddened case because higher mass galaxies tend 
to have more dust. This trend is qualitatively consis- 
tent with observations fe.g. lPapovich et al .^2001'1. It also 
leads to the counter-intuitive albeit weak trend that more 
rapidly star-forming galaxies are redder. Unfortunately, 
this correlation is not visible in the photometry (substi- 
tuting a rest-frame UV magnitude for the star formation 
rate) if the scattering effects of reddening have not been 
removed. 

- There are no obvious correlations between any pho- 
tometric measurement and the galaxies' stellar ages, 
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con sistent with the find ings of iPapovich et al.l 1)200 If) 
and iShaplev et al] 1)200 ll ^ that star formation history is 
difficult to constrain from photometry alone. 

- There is a very strong correlation between stellar mass 
and rest-frame optical magnitude, as expected. This cor- 
relation persists in the longer-wavelength IRAC bands 
(not shown) but the scatter does not decrease dramat- 
ically. We note that the correlation between optical 
flux and stellar mass in the G6 simulation at 2: = 2 
has previously been shown to agree well with observa- 
tions (,Nagamine et a l. 2005li). 

- The dashed lines in the figure of Af, vs. [3.6] show 
lines of constant mass-to-light ratio corresponding to 1 /8 
(brightest), 1/4 (middle), and 1/2 (dimmest) of the solar 
value in this band. Evidently, the stellar populations of 
B-dropouts have roughly 1/4 the solar mass-to- light ra- 
tio in rest-frame optical light. For comparison, an unred- 
dened SSP of the same metallicity with ages between 200 
and 600 Myr has rest-frame V-band mass-to-light ratios 
of 0.12 to 0.24. The difference between the two ranges 
owes entirely to reddening: As the dashed contours in 
Figure IT^ reveal, the mass-to-light ratios for the unred- 
dened galaxies fall neatly within the range given by the 
SSPs. The agreement further indicates that the galaxies' 
spectra are not dominated by bursts. 

- There is a weak trend between color and environment 
in the sense that redder galaxies are in more overdense 
areas, and the trend between stellar mass and environ- 
ment does translate into a trend between luminosity and 
environment. It will be interesting to see if stronger cor- 
relations are seen in the data; the ones predicted here 
will be quite difficult to pull out of the GOODS sample 
given additional scatter from photometric errors (which 
we have not included). 

5.3. Photometric Properties 

In this section we make predictions for observable pho- 
tometric relationships for GOODS galaxies. Validating 
these predictions would lend support to our model for 
high-redshift galaxies, while disagreements would hope- 
fully give insights into physical processes that are absent 
from the simulations. Figure 1131 shows a short-baseline 
color, a long-baseline color, and a rest-frame UV mag- 
nitude plotted against a rest-frame optical magnitude 
([3.6]). The following properties are evident: 

- There is a strong positive correlation between rest- 
frame UV and optical flux, consistent with the relation- 
ship that the more massive galaxies are more rapidly 
star- forming (Figure [TT)l . This is a fundamental predic- 
tion of our simulations, although the predicted scatter 
may be underestimated if starbursts that cannot be re- 
solved by these simulations are prevalent. 

- There is a strong correlation between the optical flux 
and color for both long-baseline (UV-optical) and short- 
baseline (UV-UV) colors. Since redder colors indicate a 
lower specific star formation rate, the simulations are re- 
pr oducing the important trend in galaxy evolution noted 
in IPapovich et al.l 1)200 If) . 

- We identify no passively-evolving, "red and dead" 
massive galaxies in our z = 4 sample. At lower red- 
shifts, such galaxies appear in a tight "red sequence" 
in color-magnitude diagrams l)Bell et al.l l20041. There is 
as yet no evidence for a red sequence at high redshifts 
0^1 although several massive, passively evolving galax- 



ies hav e been identified (Papovic h et al. 2005, in prepa- 
ratio n; lYan et alJl200l iDaddi e t al. 2005; L abbe et alJ 
120051). In any case, the color- magnitude diagrams in Fig- 
ure are clearly not bimodal. Indeed, this dearth of 
passively evolving massive systems is seen at lower red- 
shifts in models as well, and highlights an important fail- 
ure of current galaxy fo rmation models l)Somervi Ii3l200l 
iNagamine et al.l2005a|) . Therefore it is of interest to pre- 
dict where red sequence galaxies might lie on these plots 
if they existed at high redshift. 

U sing ~ 25 OOP g alaxies from the COMBO-17 Sur- 
vey, iBell et aJ] l)2004)) identified a color criterion for iso- 
lating red sequence galaxies and showed how the popula- 
tion's colors evolve from z = 0.2 1.1. One limit on the 
color of the z — A red sequence is obtained by assuming 
that the red sequence does not evol ve between z = 4 and 
z = 1.1. Using the color limit of IBell et all l)2004|) at 
z = 1.1, and converting the rest- frame U — V color and 
absolute V magnitude to the observed Veoe — [3.6] color 
and observed [3.6] magnitude, gives the upper line in the 
plot of [3.6] versus Vqo6 — [3.6]. On the other hand, since 
roughly 4 Gyr pass between z = 4 and z = 1.1 in our 
cosmology, we expect that the color of the hypothetical 
red sequence will be different at z = 4. By assuming that 
the color evolves linearly with time we obtain an upper 
limit to the amount of color evolution experienced by the 
red sequence, shown by the lower line. Both lines suc- 
ceed in isolating a sample of fairly massive red galaxies 
in the reddened sample. However, these galaxies are all 
forming stars at more than 1 Mq yr"-'^ and are red owing 
to dust, not low birthrates; if dust is removed then very 
few of the simulated galaxies are redder than either cut. 
The existence of dust-free massive red galaxies at these 
epochs, if found, would be a strong challenge for current 
models of galaxy formation. 

Figure El illustrates how ACS and IRAC data from 
GOODS may be used to constrain both the amount and 
form of reddening. We have produced photometry for 
two simulated GOODS samples from the G6 simulation, 
one in which there is no dust reddening and one in which 
there is substantial dust reddening. In the top panel 
of Figure ^1 the blue crosses represent the unreddened 
LBGs while the red contours enclose 63% and 95% of 
the reddened LBGs; in the bottom panel the red crosses 
represent the reddened LBGs while the blue contours en- 
close 63% and 95% of the unreddened LBGs. The point 
sizes are scaled to the logarithm of the stellar mass as in- 
dicated. The histograms on the side show the normalized 
distributions of colors for each sample. The unreddened 
sample consists of a dense "blue clump" along with a dif- 
fuse "red tail" of low-mass galaxies whose optical color is 
nearly independent of UV color. In the reddened sample 
the "blue clump" has been spread out along the direc- 
tion of the reddening vector — wh ich results directly fro m 
the reddening curve employed l)Calzetti et alJ 120001) — 
and many of the "red tail" galaxies have been reddened 
out of the sample. 

Roughly, the median colors can be calibrated to yield 
the typical dust extinction. Furthermore, the slope of the 
regression fit tells about the reddening law itself, since 
different laws will scatter points along different redden- 
ing vectors. As an illustration, we have included the 
regression fit to the redde ned sample for the c ase of the 
Milky Way reddening law l)Cardelli et al.lll989D . and the 
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Fig. 13. — Photometric properties of B-dropouts versus observed flux in [3.6]. (top) rest-frame UV; (second) rest-frame UV-optical color; 
(third) rest-frame UV color. The solid lines give possible color cuts for defining a high-redshift red sequence (see text). The contours 
enclose 63% and 95% of the simulated galaxies. The histograms and contours are constructed by volume-weighting the contributions of 
the D5 and G6 samples as described in the text. Rest-frame optical and UV flux show a strong correlation. 



average law for t he stars in the Small Magellanic Cloud 
(SMC) bar from 'Cord on et alJ l|2003ft . In all cases we 
have used the fiducial E{B — V) distribution to obtain 
the rest- frame V extinctions Ay- While the Calzetti 
and Milky Way laws yield similar results, the SMC law 
produces much more reddening in the rest-frame UV as 
compared to the rest-frame optical, yielding a much flat- 
ter slope in this color-color space. Although the differ- 
ences in color for the different extinction laws are some- 
what subtle, such differences may be distinguishable in 
COODS data and could help constrain the general form 
of reddening at high redshift. Unfortunately, this tech- 
nique probably would not be effective on an individual 
galaxy basis owing to photometric errors. 



5.4. High Star Formation Rate Objects 



Figure ^2 shows that the simulated COODS sample 
includes a number of galaxies with instantaneous star 
formation rates in excess of 100 Mq yr~^, including two 
in excess of ~ 1000 Mq yr~^. As noted in § 14.31 these 
galaxies are part of the bright-end excess in the rest-UV 
luminosity function (cf. Figure EJ, and may simply rep- 
resent some failing in the numerics or feedback modeling 
of high-redshift galaxies. However, the intriguing alter- 
native is that these galaxies are indeed genuine predic- 
tions of our model, and are not observed in the rest-UV 
because they are in fact much more heavily extinguished 
than our assumed dust relation (equation^ predicts. In- 
deed, such a population of heavily extinguished, rapidly 
star forming galaxies has been observed at high redshift 
as sub-millimeter galaxies (SMCs) using SCUBA. It is 
therefore of interest to investigate the physical proper- 
ties of these systems in more detail. 
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Fig. 14. — Rest-frame optical versus UV color of B-dropouts from the G6 simulation assuming fiducial reddening (red) and no reddening 
(blue). Contours enclose 63% and 95% of the galaxies. Point size scales linearly with log(M»/M0) as indicated, (top) reddened LBGs are 
given by the contour, unreddened LBGs are given by crosses; (bottom) reddened LBGs are given by crosses, unreddened LBGs are given 
by the contours. The solid, dotted, short-dashed, and long-dashed lines give regression fits to locus of B-dropouts assuming no reddening, 
fiducial reddening, a Milk Way reddening curve, and the SMC bar reddening curve. The slope of the color-color plot varies with the 
reddening curve and the median color varies with the typical amount of dust in the galaxies. 



The star formation histories of the galaxies with the 
four highest instantaneous star formation rates are shown 
in the top of Figure ^1 with the instantaneous star for- 
mation rate given in the top left corner and the stellar 
masses given above each plot. Despite their prodigious 
star formation rates, they are not undergoing significant 
starbursts; their star formation rates are all increased 
by a modest 2 — 3x over the rate given by equation [S| 
Closer inspection reveals that all of these galaxies reside 
in highly overdense regions resembling the cores of col- 
lapsing clusters, and owing to the large number of smaller 
galaxies in the vicinity, they exhibit a complex morphol- 
ogy that could be construed as indicating interaction. 
However, none of them appear to be the merger of two 
massive gas-rich galaxies. While our simulations cannot 
resolve the detailed gas dynamics within each galaxy, it 



seems unlikely that their star formation is being driven 
by dynamical instabilities that funnel gas to a small cen- 
tral core as in merger simulations (e.g. Barnes & Hern- 
quist 1991, 1996). Instead, it appears likely that star 
formation is occurring throughout the galaxy, and the 
high SFRs are purely a result of the high gas accretion 
rates expected for these massive systems. 

For more quantitative insight, we compared their in- 
stantaneous SFRs with their 100 Myr-averaged SFRs 
(bottom left) and their birthrates (bottom right). The 
two most prodigious star formers are forming more than 
1000 M0 yr~^, or about three times their 100 Myr av- 
erages. They show birthrates in excess of 5 (i.e., they 
are forming stars over 5 times more rapidly than would 
be needed to acquire their stellar mass at z = 4 within 
a Hubble time). The other two galaxies are forming 
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Fig. 15. — (top) Star formation histories of the four most rapidly star forming galaxies in the simulated samples. The instantaneous 
star formation rates and stellar masses at 2 = 4 are given in the top left; the star formation histories are binned in intervals of 10 Myr. 
(bottom left) Instantaneous versus time-averaged star formation rate for the GOODS-selected galaxies in the G6 simulation, with the four 
most rapid star formers marked by triangles. The straight lines give the location of galaxies whose instantaneous star formation rates are 
(bottom line) equal to and (top line) twice the time-averaged rates, (bottom right) Instantaneous star formation rate versus birthrate for 
the same sources. The lower envelope arises from our 10^'^^Mq mass resolution cut. These massive, rapid star formers are experiencing 
at most a mild burst of star formation and may be identifiable with high-redshift submillimctcr galaxies. 



stars at 500-600 Mq yr~^, or roughly twice their 100 
Myr averages, and show instantaneous birthrates of w 
3-4. In fact, the instantaneous versus averaged star for- 
mation rate plot suggests that the galaxies whose SFRs 
exceed 100 Mq yr~^ depart smoothly from the trend set 
by the less rapidly star-forming (M^/M© < 100) galax- 
ies. Hence these high-SFR galaxies are not anomalies, 
but a natural extension of the trends from lower masses 
and SFRs. 

The regime of > 1000 Mq yr-i objects has gener- 
ally been thought to be populated by SCUBA sources, 
which are in turn ultraluminous IR sources thought 
to be scaled-up analogs of today's ULIRGs. Can our 
high-SFR galaxies be SCUBA sources? Unfortunately, 
very little information is available on sub-mm galaxies 



at 2 ~ 4. However, we can compare to samples at 
z ^ 2 — 3, and the results are intriguing. First, the 
space density of our two SFR> 1000 Mq yr~^ galax- 
ies is comparable to that of z ~ 2.5 sub-mm galaxies 
with Lfir > IO^/i-^Lq ijChaoman et al.ll2005j) . Hence 
they have broadly similar star formation rates and num- 
ber densities, even if the numbers of SMGs decreases 
by a factor of a few from z ~ 2.5 4 (as seen in 
Figure 4 of Cha pman ct al. 20Q3)- Even more curious 
is that IChapman et alJ lj200^ find, for a sample of 12 
sub-mm galaxies with z « 2.2, the majority show ex- 
tended radio emission indicative of a spatially-extended 
starburst. This seems more in accord with the nature of 
our simulated objects than local ULIRGs that look more 
like canonical major mergers with bright central star- 
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bursts, although we note that the merger-driven model 
of Hopkins et al. (2005a, b) for the simultaneous origin of 
starbursts and quasars also predicts a density of sub-mm 
galaxies similar to that observed (Hopkins et al. 2005c). 

It would clearly be interesting to compare the statis- 
tical properties of the simulated massive galaxies to ob- 
servations in more detail. Although th is is beyond the 
scope o f the present work, we note that iNagamine et al.l 
(j2005b|) . who used the same SPH simulations, found good 
agreement between the simulated and observed number 
densities of massive galaxies at z =1-3 and used this 
finding to argue that hydrodynamic simulations do not 
suffer from the "mass-scale problem." Since our massive 
galaxies at z — A must be the progenitors of their sample, 
their work lends further indirect support to the sugges- 
tion that the our massive galaxies may correspond to an 
observable population such as the SMGs. 

In summary, our simulations produce several rapidly 
star-forming objects whose properties are broadly in ac- 
cord with those of SMGs seen with SCUBA. They are 
forming stars at a rate that is at most mildly in excess 
of what is expected for their stellar mass, indicating that 
gas infall and not merger-induced dynamical instabilities 
is driving their star formation. For these objects to be 
SMGs, their dust extinction must lie significantly above 
the relation we assume for the bulk of our LEG popu- 
lation (equation^, but given the ad hoc nature of that 
relation, this does not seem far-fetched. More detailed 
model ing is required, such as that done in lFardal et al.l 
to see if this hypothesis is viable. 

6. CONCLUSION 

We have analyzed the physical and photometric prop- 
erties of galaxies at z = 4 in a cosmological hydro- 
dynamic simulation of a ACDM universe. We employ 
a novel prescription for modeling dust reddening based 
on each simulated galaxy's metallicity, and determine a 
mass criterion that robustly resolves the star formation 
history of these galaxies. We examine the properties of 
galaxies that satisfy the selection criteria for the z ^ 4 
B-dropout sample of the GOODS survey. 

We subdivide our conclusions into two broad areas: 
Observed properties of the GOODS sample, and physical 
properties of high-redshift galaxies. For the former, we 
find: 

• The simulated luminosity density from rest-frame 
UV through rest-frame optical bands is in broad 
agreement with observations. While the normal- 
ization is uncertain to ~ 2 x owing to dust redden- 
ing, the shape of the luminosity density curve is 
insensitive to reddening, for reasonable dust mod- 
els. We find that the simulated luminosity density 
curve includes a conspicuous 4000 A break that is 
not seen in observations; whether or not this dis- 
crepancy is genuine can be tested with increased 
areal coverage in the iJ-band, as well as luminosity 
density measurements in IRAC bands. 

• The simulated rest-frame UV luminosity function is 
also in broad agreement with the data. The 
turnover at the faint end owes entirely to the mag- 
nitude limit; the intrinsic galaxy population in the 
simulation has a steep faint-end slope of a ^ —2. 



The simulated rest-frame optical luminosity func- 
tion is similar to the rest-UV luminosity function, 
shifted brighter by 1-2 magnitudes. 

• The simulated rest-UV LF shows an excess of 
galaxies with 1775 < 25 that could indicate a failing 
of the simulation, but may also represent a popula- 
tion of highly-obscured, rapidly star- forming galax- 
ies that are not modeled correctly by our reddening 
prescription. 

• The GOODS B-dropout selection criteria miss low- 
mass galaxies owing to the brightness cut at zgso 
and some massive galaxies owing to the blueness 
cut in Vgoe ~ -2^850- Assuming our fiducial reddening 
prescription, we find that B-dropout samples ob- 
serve less than 50% of the stellar mass density in 
the univer se at this epoch. Th e RJL color criterion 
defined bv lDaddi et al.l l|2004D does not exclude any 
of our simulated galaxies for redshifts below 4.7. 

For the physical properties of simulated high-redshift 
galaxies, we find: 

• Stellar mass and star formation rate are tightly cor- 
related (equation |SJ), leading to the fundamental 
prediction that GOODS B-dropouts are the most 
massive galaxies to have formed by z ^ 4. 

• The B-dropout stellar masses fall predominantly in 
the range log(M*/M0) = 9-10 with a significant 
tail of massive galaxies out to log(M*/M0) « 11, 
in good agreement with recent observational assess- 
ments and somew hat more massive than the semi- 
analytic sample of lldzi et alJ l|2004D . 

• B-dropouts are metal-enriched and obey a tight 
stellar mass-metallicity relation. The slope of this 
relation is in excellent agreement with the low- 
redshift relation of Trcmonti et al. ( 2004) and the 
normalization is such that B-dropouts of a given 
stellar mass are roughly 0.6 dex less enriched than 
a typical low-redshift galaxy of the same mass. 

• The median ages of the star particles in the B- 
dropouts range from 200-600 Myr, consistent with 
observational constraints. However, stellar age is 
not tightly correlated with any other physical or 
photometric parameter. 

• There is a weak correlation between stellar mass 
and environment in the sense that the most mas- 
sive galaxies are found in the densest environments; 
more careful modeling is required to determine if 
this is observable in the GOODS data. 

• A significant population is forming stars at rates 
in excess of 100 Mq yr~^, including several exceed- 
ing 1000 Mq yr~^. These high star formation rates 
are not dramatic excursions over a slower quiescent 
rate, but rather are simply the rates that are ap- 
propriate for their large stellar masses. The most 
rapid star formers have properties intriguingly sim- 
ilar to SCUBA galaxies, though a firm assessment 
of this hypothesis will require further work. 
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• Rest-frame UV and rest-frame optical fluxes are 
both correlated with star formation rate; however, 
owing to the tight correlation between stellar mass 
and star formation rate together with the scatter- 
ing effect of dust reddening on rest-frame UV flux, 
rest-optical fluxes are more tightly correlated with 
star formation rate than rest-UV fluxes, even if 
dust reddening can be accurately removed. 

• More massive galaxies are redder. This trend is 
present even in the absence of dust reddening, al- 
though dust enhances it since more massive galax- 
ies are dustier in our reddening prescription. This 
results in rest-frame optical flux being correlated 
with color in the sense that brighter galaxies are 
redder, consistent with recent observations. 

• Our simulations do not produce a bimodal distribu- 
tion of galaxy colors as seen in the red sequence at 
low redshift. We have predicted \ Yhere the d i vision 
might occur if the trend noted bv lBell et all 1)20041) 
at z < 1 persists out to 2 ~ 4 and find that many 
dusty star-forming galaxies may be red enough to 
satisfy a red sequence color cut, but this sequence 
would not be tight. 

• The distribution of galaxies on a UV-to-optical 
color-color plot may be used to constrain the typ- 
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ical amount and form of reddening sulfered by B- 
dropouts although it probably cannot be used on 
individual galaxies. 

In summary, we find that current simulations of galaxy 
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observations at z ~ 4, and provide valuable insights into 
the nature of high redshift galaxies. We look forward to 
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to test galaxy formation models in more detail. 
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